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FOREWORD 


As it has in centuries past, the Earth's magnetic field still plays 
a vital role in global navigation. All navigational aids or 
attitude/heading reference systems (AHRS) , regardless of their 
operating principles, must speak a common language. That common 
language is in terms of the Earth's magnetic declination. 
Consequently, magnetic-related navigational aids are integrated, 
in the form of computer hardware and software, into virtually every 
major weapons system of the Army, Air Force, Navy, and Marines. 
In order to maintain optimum performance, these systems must be 
periodically updated with regard to the Earth's magnetic field, 
which is a dynamic entity that changes slowly but erratically with 
time . 

For well over a century, it has been the responsibility of the U.S. 
Naval Oceanographic Office to monitor the Earth's changing magnetic 
field and periodically report on these changes in the form of 
magnetic charts and mathematical models. For the past forty years, 
this task has involved an intensive data collection effort through 
the Project MAGNET program, which in April 1990 made the transition 
from primarily aeromagnetic surveying to satellite surveying with 
the launch of the Polar Orbiting Geomagnetic Survey (POGS) satellite. 
Follow-on satellite missions to secure data for future needs, well 
into the twenty-first century, are now being vigorously pursued. 

This report is a comprehensive summary of the cooperative effort 
between the U.S. Naval Oceanographic Office and the British 
Geological Survey in producing the 1990 Epoch World Magnetic Model, 
WMM-90. 
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SECTION 1. THE GEOMAG ALGORITHM AND THE 1990 MODEL 
1 . 0 Introduction 

The Earth's magnetic field, as measured by a magnetic sensor on or 
above the Earth's surface, is actually a composite of several magnetic 
fields generated by a variety of sources, which are superimposed on 
each other and which interact with each other. The most important 
of these geomagnetic sources are: 

a. the Earth's fluid outer core; 

b. the Earth's crust/upper mantle; 

c. the ionosphere; and 

d. the magnetosphere. 

The magnetic variation algorithm (GEOMAG) is a Fortran subroutine 
which is based on a spherical harmonic expansion of the Earth s 
magnetic field, the coefficients of which comprise the World Magnetic 
Model (WMM) . These coefficients are produced jointly by the U.S. 
Naval Oceanographic Office (NAVOCEANO) ' s Geopotential Division and 
the British Geological Survey (BGS) . The WMM is distributed by 
NAVOCEANO for the Defense Mapping Agency (DMA) in accordance with DMA 
Instructions 8000.1 and 8000.2. The WMMs are usually produced at 
5-year intervals and are composed of two parts: a main field model, 

which describes the Earth's magnetic field at some base epoch, and 
a secular variation model, which accounts for the slow temporal 
variations in the main geomagnetic field from the base epoch to a 
maximum of 5 years beyond the base epoch. For example, the base epoch 
of the WMM-90 magnetic field model is 1990.0. This model is therefore 
considered valid between 1990.0 and 1995.0 and will subsequently be 
replaced at 1995.0 by the WMM-95 magnetic field model. 

It is extremely important to recognize that the WMM series of 
geomagnetic models and the charts produced from these models 
characterize only that portion of the Earth's magnetic field whi^h 
is generated by the Earth's fluid outer core. The portions of the 
geomagnetic field generated by the Earth's crust, mantle, ionosphere, 
and magnetosphere are not represented in these models. Consequently, 
a magnetic sensor such as a compass or magnetometer may observe spatial 
and temporal magnetic anomalies when referenced to the appropriate 
WMM. In particular, certain local, regional, and temporal magnetic 
declination anomalies can exceed 10 degrees. Anomalies of this 

magnitude are not common, but they do exist. Declination anomalies 
on the order of 3 or 4 degrees are not uncommon, but are of small 
spatial extent and relatively isolated. On land, spatial anomalies 
are produced by mountain ranges; ore deposits; ground which has been 
struck by lightning; geological faults; and cultural features such 
as trains, planes, tanks, railroad tracks, power lines, etc. In ocean 
areas, spatial anomalies are produced by continental margins, 
seamounts, oceanic ridges, trenches and fault zones, and ships and 
submarines. Temporal anomalies in either ocean or land areas can 
last from a few minutes to several days and are produced by ionospheric 
and magnetospher ic processes which are driven by the solar wind. 



Magnetic storms in particular can cause severe and persistent magnetic 
anomalies. Even in periods of quiet solar activity, significant 
spatial and temporal magnetic anomalies are found in the polar and 
equatorial regions of the Earth, where magnetic fields produced by 
ionospheric current systems, such as the auroral electro jets and the 
equatorial electro jet, are always present. Most of the possible 
sources of magnetic anomalies are comparatively isolated in either 
space or time. Therefore, from a global perspective, the 
root-mean-square (RMS), declination (DEC), and inclination (DIP) 
errors at sea level of the WMM are estimated to be less than 0.5 
degrees in ocean areas and less than 1.0 degrees in land areas at the 
Earth's surface over the entire 5-year life of a particular model. 
Also, the RMS errors at sea level of the horizontal (H) and total 
intensity (F) components of the WMM over ocean and land areas are 
estimated to be less than 200 nanoteslas (nT) over the entire 5-year 
life of the models. 


1 . 1 The Math ematical Model 

The Earth's magnetic field has associated with it a geomagnetic 
potential V'(r,e,<}),f), which can be expressed in spherical coordinates 
in terms of a spherical harmonic expansion of the following form: 


V(r,e,<M) 



Re 

r 


Y + i 

/ 


I {£«n( x ) cosm< l> + sin/n<|>}7 > "(9) 

m *0 


(i) 


where the spherical coordinates <r, correspond to the radius from 
the center of the Earth, the colatitude (i.e., 90* - latitude), and 
the longitude. R E is the mean radius of the Earth; g m (x) and h m {x) are 
referred to as the Gauss coefficients at time X, where X is the time 
in years (e.g., 1987.312). P?(9) represents a particular associated 

Legendre polynomial of degree n and order m. These polynomials are 
functions of the colatitude 0. The Gauss coefficients are slowly 
varying functions of time and are expressed in the form: 

£/vn( X ) = 8 mSJ EPOCH) + £ wn( X_ ^ EPOCh) (2 a) 

^un(t) = h wJJ epoch) + h^ix — T epoch) (2b) 


where Tepoch is the base epoch of the model, which for WMM-90 is 1990.0. 
Thus, 8mSTepoch) an< 3 hnnStT epoch) ar e the Gauss coefficients of the WMM at 
the model's base epoch, while and h m (pronounced g m dot and 
dot) are the annual rates of change of the Gauss coefficients. The 
Gauss coefficients g^EPocH) and h m (J epoch) and their annual rates of 
change are spherical harmonic coefficients. The Gauss coefficients 
gim(T epoch) a ^d H^Tepoch) characterize the Earth's main magnetic field at 
the base epoch of the model, Tepoch/ while and h ^ characterize the 
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secular change of the- Earth's main magnetic field during the 5 year 
life of the model. These coefficients, up to degree and order 12 for 
the main field and up to degree and order 8 for the secular variation 
of the main field, comprise the WMM . Currently, the secular variation 
model from degree 8 through degree 12 is padded with zeros. 

The Earth's magnetic field B(r, 0,<j>,T) is a vector quantity having three 
components which correspond to the projection of the magnetic field 
vector onto the three coordinate axes. Thus, 5 r (r,0,<t>,t) is that portion 
of the field pointing in the radial direction (i.e., perpendicular 
to the surface of the Earth), B 8 (r,0,<j>,x) is that portion of the field 
pointing locally due south, and B^(r,Q,^,x) is that portion of the field 
pointing locally due east. The magnetic field vector can be computed 
from the geomagnetic potential by taking its gradient, thus. 

fi(r,e,<M) = -W(r,0,$,T) 


Consequently, the magnetic field components are related to the 
geomagnetic potential as follows: 


B r O\0,<J>,t) = 
S e (r,0,< M) = 
^(r,0,<M) = 

which yield the 
B r (r,0,<J>,T) = 

S e (r,0,<M) = 
^(r,0,<M) = 


aV(r,0,( M) 

dr 

1 0V'(r,0,4>,x) 

r 00 


1 0 V (r , 0, <|>, t) 


r sin0 


0<J> 


following spherical harmonic expansions: 


X (n + 1) 

H * 1 


/fl r V 2 - 


V r ) 


X {^(x) cos m<t> + h^x) sin m4>}P"(0) 


N 

- X 

n = 1 


/» \» + 2 


dP?(Q) 

y I X o {^ WB (x) cos m<j> + h m (x) sin m <J> } — ^ 


1 N 

—— X 

sin 0 » = i 


^ r ) 


X m{g )WI (x)sin/«<|>-/i M ,(x)cos/n<!)}/ , "(0) 


m = 0 


(4a) 

(4b) 

(4c) 

(5a) 

(5b) 

(5c) 


It must be noted that the Gauss coefficients g m (x) and h^x), as well 
as the associated Legendre polynomials and their derivatives, are 
Schmidt normalized by an international agreement (circa 1930) of the 
International Union of Geodesy and Geophysics. This particular 
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normalization allows one to determine which terms of the spherical 
harmonic model are the most significant simply by a cursory inspection 
of the model coefficients. The Schmidt-normalized associated Legendre 

Polynomials P"(0) are related to the unnormalized associated Legendre 

Polynomials P' u "(0) (note position of indices) by the following relation: 


p;(0) = 


The Schmidt normalization factors 5*" and the unnormalized associated 
Legendre Polynomials P""( 0) are computed via recurrence relationships 
as follows: 

P°°(Q) = 1 (7a) 


p“(0) = 

sin 0P"' 1,m_1 (0) 


m =n *0 

(7b) 

P' un (Q) = 

COS0P'" 1 '" - k W, P"' 2 '"’ 


m * n,n > 1 

(7c) 

dP°°(9) 

d9 

0 



( 7d) 

dP m (0) 
d9 

sin 9 dp -;:" (e) + cos9/>” 
at) 


,m=n* 0 

(7e) 

dP*”(9) 

dd 

c°s0 — — — - sin0P" _1 "(0) 
at) 

l/U ndP"- 2 ’ m 

de 

,m *n,n > 1 

( 7 f ) 


where : 

= (n - l) 2 -m 2 (8) 

(2n - l)(2n -3) 

and where it is understood that the undefined polynomials P" 1,o (0) and 
— — -(0) are to be set equal to zero. Similarly, 

ufj 

£<x> j (9a) 

S” = ±2 — i 5*- 1 - 0 >0 ( 9 b ) 
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(9c) 


S 


nm 



(n-m + iy 
n+m 


S 


n,m- 1 


{J = 2 for m = lj 
j/ = l for m> 1J 


Also, computed via recursion relations are the longitudinally dependent 
functions cos (m§) and sin(m4>), which are computed as follows: 


sin(m <+)) = 0 ,m—0 

cos(m<))) = 1 ,m =0 

sin(m4>) = sin(<j>)cos(m - 1)<|) + cos((|)) sin(m - 1)<> 

cos(m<))) = cos(<J))cos(m - 1)<J> - sin(4>)sin(m - 1)<> 

1.2 Coordinate Transformations 


,m> 0 
,m > 0 


(10a) 

(lOt) 

(10c) 

(lOd) 


GEOMAG is intended to compute various components of the geomagnetic 
field in a geodetic coordinate system that uses the WGS-84 ellipsoid 
as the reference ellipsoid. However, the mathematical analysis in 
the previous section is based on spherical coordinates. Consequently, 
some coordinate transformations are necessary. A three-step procedure 
is required. 


a. Convert the geodetic latitude, longitude, and altitude 
to spherical coordinates ( r, 9, <J>) . 


b. Compute the magnetic field components B,(r,0,({),T), fl 9 (r,0,<(),x), and 


c. Rotate the magnetic field components from spherical 
coordinates to geodetic coordinates yielding the magnetic field 
components B x (k,§,h,X), B Y M,h,x) f andB z M,h,x), which are the projections 

of the magnetic field vector B(k,ty,h,x) onto the X-north, Y-east, and 
Z-vertically down coordinates of a local rectangular coordinate system 
defined by the tangent plane to the ellipsoid which is concentric 
about the WGS-84 reference ellipsoid but which encompasses the point 

The transformations in step a are as follows: 


COS0 = 

sinX 

(11a) 

V Q 2 cos 2 X + sin 2 X 

sin0 = 

V 1 - cos 2 0 

(lib) 

where, if a 
of the WGS- 

and b are respectively the 
8 4 ellipsoid: 

semi-major and semi -minor axes 


s 



Q 


( 12 ) 


Wa 2 -(fl 2 -ft 2 )sin 2 X + a 2 
Wa 2 -(a 2 -6 2 )sin 2 X + £ 2 

Furthermore : 


r 


2 


/i 2 + 2/Wa 2 (a 2 -£ 2 )sin 2 X 


g 4 -(a 4 -fr 4 ) sin 2 K 
a 2 -(a 2 -& 2 )sin 2 X 


(13) 


The transformation in step c depends on the rotation angle a through 
which the magnetic field vector must be rotated in going from spherical 
to geodetic coordinates. This rotation angle is defined by the 
following rotations: 

cos a = { h + Va 2 cos 2 X + Z? 2 sin 2 X, }/r 


sin a = (a 2 — b 2 ) cos^sinX/ {r i/a 2 cos 2 X + 6 2 sin 2 X } 


a = 



(14c) 


Consequently, the components of the magnetic field vector in geodetic 
coordinates may be computed as follows: 


B x (k,$,h,x) = - cosaS e (r, 9,$,x) - sin a B r (r ,0,$,x) 


B z (k,$,h,x) = sin a 5 e (r , 0, <(>, x) - cos a B r (r , 0, <|», t) 


(15a) 

(15b) 

(15c) 


From these rectangular components of the geomagnetic field, it is 
possible to construct all others. In particular, the following 
parameters may be computed: 


B H (k,ty,h,X) 

B F (Kty,h,x) 

BpChj h , x) 


^B$(K$,h,x)+B}M,h,x) 

(Horizontal 

Intensity) 

(16a) 


(Total 

Intensity) 

(16b) 

-i ByQ^y ^ h , x) 1 

tan 1 1 : — rr 

(Declination) 

(16c) 

{B x M,h,x)\ 
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(Inclination) 


( 1 6 d) 




, \ B z M,h,x) \ 


B c (k,§,h,x) 


[B d -* 

K+4> 


x>ol 

x<o] 


(Grid 


Variation) 


<16e) 


1 .3 The Computer Al gorithm 

The Gauss coefficients at the base epoch, T'epoch/ 3re stored in array 
C so that the lower half of array C is occupied by the even harmonic 
Gauss coefficients g m (J epoch) t while the upper half of array C is 
occupied by the odd harmonic Gauss coefficients B^JJ epoch) • Table 1 
illustrates the details of the storage scheme, which is equivalent 
to the following mathematical assignments: 


8 

h 


nm 

m,n + 1 


,m <n 
,m> n 


which implies that: 

8 run = ,m <, H 

h m = C„_ M ,m<,n,m* 0 


(17) 


(18a) 

(18b) 


The annual rates of change of the Gauss coefficients are stored in 
array CD (which stands for C) so that the lower half of array CD is 
occupied by the even harmonic coefficients while the upper half 
of the array is occupied by the odd harmonic coefficients /i* n - Table 
2 illustrates the details of the storage scheme for array CD . It is 
essentially the same as table 1 for array C and corresponds to the 
following mathematical assignments: 


§ nm 

B m, n + l 


which implies that: 

& nm ^ nm 

B nm ~ ^m-l.n 


,m <i n 
,m> n 


,m <>n 

,m <n,m * 0 


(19) 


(20a) 

(20b) 
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TABLE 1. ARRANGEMENT OF MAIN FIELD COEFFICIENTS IN ARRAY C m 


n\m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0 

£oo 

*11 

^21 

/l 3l 

/z 4 i 

HI 

*61 

*71 

^81 

^91 

^10,1 

^n,i 

^12,1 

1 

8 io 

£n 

ra 

*32 

IS 

ra 

*62 


^82 



^92 

^10,2 

^11,2 

^12,2 

2 

8 20 

£21 

8 22 

IS 

^43 

^53 

*63 

*73 


*93 

^10,3 

^11,3 

^12,3 

3 

8 30 

£31 

8 32 

£ 33 

A44 

^54 

*64 

*74 


^94 

ra 

1,4 

^12,4 

4 

8 40 

£41 

£42 

£43 

£44 

*55 

h 6 s 

*75 

^85 

^95 

ra 

^11,5 

^12,5 

5 

#50 

£51 

£52 

£53 

£54 

£55 

*66 

*76 

^86 

^96 

ra 

^11,6 

^12,6 

6 

<?60 

<?61 

£62 

£63 

£64 

£65 

£66 

*77 


h<n 

^10,7 

^11,7 

^12,7 

7 

#70 

£71 

£ 72 

£73 

£74 

£75 

£76 

£77 

□ 

^98 

^10,8 

^11,8 

^12,8 

8 

£80 

<?81 

£82 

£83 

£84 

£85 

£ 86 

£#7 

£88 

*99 

^10,9 

^11,9 

^12,9 

9 

£ 90 

£91 

£ 92 

£93 

£94 

£95 

£96 

£97 

£98 

£99 

ra 

^11,10 

^12,10 

10 

£io,o 

810,1 

£10,2 

£l0,3 

£10,4 

£l0,5 

£l0,6 

£l0,7 

£l0,8 

£10.9 

£10, 10 

*11,11 

^12,11 

11 

£n,o 

£11,1 

£11,2 

£ll,3 

£11,4 

£ll,5 

£u,6 

£ll,7 

£ll,8 

£11,9 

£11,10 

£ll,ll 

^12,12 

12 

£l2,0 

£12,1 

£12,2 

£12,3 

£l2,4 

£l2,5 

£12,6 

£12,7 

£ 12,8 

£ 12,9 

£12,10 

£l2,ll 

£l2,12 





























































































































































TABLE 2. ARRANGEMENT OF SECULAR VARIATION COEFFICIENTS IN ARRAY C m 


n\m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0 

§ oo 

*11 

^21 

^31 

h 4l 

*51 

*61 

/z 7 i 

^81 

fl 9x 

h 10. 1 

*n,i 

^ 12, 1 

1 

§10 

§n 

KB 

^32 

m 

^52 

IB 

IB 

^82 

^92 

^ 10,2 


^ 12,2 

2 

& 20 

5 21 

§22 

^33 

IB 
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The numerical values of the Gauss coefficients at the base epoch and 
their corresponding annual rates of change for the WMM-90 geomagnetic 
model are listed in table 3. These numerical values are inserted 
into arrays C and CD through data statements. The base epoch of the 
model is also assigned through a data statement. In order to update 
the GEOMAG algorithm to a new epoch geomagnetic model such as WMM-95, 
it is necessary to replace only the data statements with the new model 
coefficients and the new base epoch. 

Important parameters in the GEOMAG routine and their mathematical 
correspondences are: 

A ~a = 6378.137 km 
B~b = 6356.7523142 km 
RE~R e = 6371.2 km 
TIME ~ x 
EPOCH ~ T epoch 

DT ~ x - T epoch 
ALT ~ h 

SNORM(N,M ) ~ 5"* 

K(N,M) ~ k" m 
GLAT - X 
GLON ~ <t> 

SP(M) ~ sin(m<t>) 

CP{M) ~ cos (m<^) 

ST ~ sin(9) 

CT ~ cos(0) 

CA ~ cos(a) 

SA ~ sin(a) 

BR ~ B r 
BT ~ B d 
BP ~ B+ 

BX ~ B x 
BY ~ By 
BZ ~ B z 
BH ~ B h 
DEC ~ B d 
DIP ~ B, 

TI ~ B ti 
MAXDEG ~ N 
MAXORD ~ M =N 

P{N,M ) ~ />*" 

DP(N,M) ~ dJ £ 
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TC ~ C + (t— Tzpqch) C 
CD ~ C 
Q2 ~ Q 2 


Note that R E is not intended to be the mean radius of the WGS 84 
ellipsoid. It is the mean radius of a modified IAU-66 ellipsoid. 


The GEOMAG algorithm is organized into two modules, each with its own 
entry point. The first is an Initialization Module. Its purpose : a 
to compute all constants such as the recursion relation factors for 

the associated Legendre polynomials k" , the Schmidt normalization 

factors S"", and any other parameters that do not depend on position 
or time. The entry point for this module is GEOMAG (MAXDEG) . The 
parameter MAXDEG determines the maximum degree and order of the 
magnetic model to be used in the computations. Normally, MAXDEG - 
12, which is the maximum degree and order of the WMM series geomagnetic 
models In order to reduce computation time, MAXDEG may be set to 
a number less than 12 (e.g., 8 or 10). However, the accuracy of the 
computed magnetic parameters is correspondingly reduced. MAXDEv must 
be set in the calling program. The second module is the Processing 
Module, which has the entry point 


GEOMG1 (ALT, GLAT, GLON, TIME, DEC, DIP, TI, GV) . 


The purpose of this module is to compute the magnetic decimation, 
inclination, total intensity, and grid variation of each geodetic 
position and time supplied to it. The units of the parameters in the 
argument list of the GEOMG1 entry point are as follows: 


ALT ~ kilometers (e.g., 5.314) (In) 

GLAT ~ degrees (e.g., 33.716) (In) 

GLON ~ degrees (e.g., -163.315) (In) 

TIME ~ years (e.g., 1992.427) (In) 

DEC ~ degrees (e.g., -121.734) (Out) 

DIP ~ degrees (e.g., 48.387) (Out) 

TI ~ nanoteslas (e.g., 35781.7) (Out) 

GV ~ degrees (e.g., 51.768) (Out) 


The computed magnetic field parameters are referenced to the WGS-84 
ellipsoid. The last parameter, GV, is the grid variation which is 
computed only in the polar regions (i.e., above + 55* latitude or 
below - 55' latitude). Outside of this region, a value of 999.0 is 
dummied in. It is referenced to grid north of a polar stereographic 
projection. The model is considered valid at altitudes ranging from 
sea level to 1000 km. 
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TABLE 3 . WMM -90 SCHMIDT NORMALIZED GAUSS COEFFICIENTS ( nT ) 


n 

m 

8 urn 

Km 

8 nm 

h*m 


1 


0 

1 

0 

1 

2 

0 

1 

2 

3 

0 

1 

2 

3 

4 
0 
1 
2 

3 

4 

5 
0 
1 
2 

3 

4 

5 

6 
0 
1 
2 

3 

4 

5 

6 
7 
0 
1 
2 

3 

4 

5 

6 

7 

8 


- 29780.5 

.0 

16 . 0 

.0 

- 1851.7 

5407.2 

9 . 3 

- 13.8 

- 2134.3 

.0 

- 11.7 

.0 

3062.2 

- 2278.3 

3.7 

- 12.8 

1691.9 

- 384.3 

1.8 

- 14 . 9 

1312.9 

.0 

2 . 1 

.0 

- 2244.7 

- 284.9 

-7 . 6 

3 . 1 

1246.8 

291.7 

.0 

.8 

808 . 6 

- 352.4 

- 5.8 

- 11.3 

933.5 

.0 

-.8 

. 0 

784 . 9 

249.4 

1 . 0 

3.3 

323.5 

- 232.7 

- 7.4 

3.7 

- 421.7 

91.3 

.8 

2.8 

139.2 

- 296.5 

- 6.4 

.0 

- 208.3 

.0 

1.7 

.0 

352.2 

40.8 

.0 

.0 

246.5 

148.7 

.0 

-2 . 1 

- 110 . 8 

- 154.6 

- 2.7 

1.2 

- 162.3 

- 67.6 

.0 

1.2 

- 37.2 

97.4 

3.0 

. 6 

59.0 

.0 

.8 

.0 

63.7 

- 14.7 

.0 

-.6 

60.0 

82.2 

1.5 

-.6 

- 181.3 

70.0 

.0 

.0 

.4 

- 56.2 

.0 

- 2 . 3 

15.4 

- 1.4 

.0 

.0 

- 96.0 

24.6 

.0 

.0 

76.1 

.0 

.5 

.0 

- 62.1 

- 78.6 

.0 

.6 

1.3 

- 26.7 

-.9 

.8 

30.2 

.1 

1.5 

.0 

4.7 

19.9 

2.7 

.0 

7.9 

17.9 

- 1.0 

.0 

10 . 1 

- 21.5 

.0 

.4 

1 . 9 

- 6.8 

.0 

.0 

22 . 9 

.0 

.0 

.0 

2.3 

9.7 

-1 . 1 

. 4 

- 1.2 

- 19.3 

.0 

- . 8 

- 11.7 

6.6 

.0 

. 5 

- 17.5 

- 20.1 

-2 . 1 

. 3 

2.2 

13.4 

.0 

. 5 

5.7 

9.8 

1.0 

.0 

3.0 

- 19.0 

.0 

-.7 

- 7.0 

- 9.1 

.0 

.0 
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SECTION 2. THE 1990 EPOCH WORLD MAGNETIC MODEL (DERIVATION) 
2 . 0 Overview 


There were four major data sets available for the 1990 model. These 
were: the MAGSAT satellite data collected during 1979 and 1980; the 
DE-2 satellite data collected from 1981 through 1983; Project MAGNET 
aeromagnetic data collected between 1980 and 1990; and geomagnetic 
observatory annual magnetic means data collected between 1980 and 
1990. The global distribution of these data is illustrated in charts 
1 through 4 . 

Four factors which affect the quality of the model produced and which 
influence the overall approach taken to produce the model are: 

a. The age of the data relative to the model epoch; 

b. The temporal coherence of the data; 

c. The spatial uniformity of the data; and 

d. The data density. 

With respect to these factors, none of the four data sets are ideal. 
All four data sets, especially the satellite data sets, are dominated 
by older data. The Project MAGNET data, in addition, are neither 
temporally coherent nor spatially uniform. Furthermore, the 
observatory annual means data are sparse and suffer from severe spatial 
nonuniformity . 

The modeling objective is to create two spherical harmonic models. 
One model characterizes the Earth's main (core-generated) magnetic 
field at the 1990.0 epoch. The other model characterizes the Earth's 
secular (slow temporal) magnetic variations of Earth core origin for 
five years beyond the 1990 epoch. 

Given the objective and the available data, the following procedure 
was adopted: 

a. Use the observatory annual magnetic means to create two 
definitive secular variation models, the first covering the 5-year 
interval 1980 to 1985, and the second covering the 5-year interval 
1985 to 1990 . These are referred to as the 1982.5 and 1987.5 definitive 
secular variation models, respectively. 

b. Use the observatory annual magnetic means to create, by 
extrapolation, one predictive secular variation model covering the 
5-year interval 1990 to 1995. It is referred to as the 1992.5 
predictive secular variation model. 

c. Use the two definitive secular variation models to push the 
satellite and aircraft magnetic field observations forward or backward, 
as appropriate, to 1985.0. 
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CHART 1 . PROJECT MAGNET DATA DISTRIBUTION (FROM SURVEYS 
PERFORMED DURING THE PERIOD 1980-1989) 
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CHART 2. MAGSAT DATA DISTRIBUTION 
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CHART 4. GEOMAGNETIC OBSERVATORY DISTRIBUTION 






d. Create a 1985.0 epoch main field model using the time-ad jusced 
satellite and aircraft magnetic field observations via a weighted 
least-square inversion. 

e. Use the 1987.5 definitive secular variation model to push the 
spherical harmonic coefficients of the 1985.0 epoch main field model 
forward to the 1990.0 epoch, thereby yielding the 1990.0 epoch main 
field model . 

f . Combine the 1990 epoch main field model coefficients with the 
1992.5 predictive secular variation model coefficients to form the 
1990 World Magnetic Model, WMM-90. 

A by-product of this procedure is a revised 1985.0 epoch World Chart 
Model which is obtained by combining the 1985.0 main field coefficients 
generated in step d with the 1987.5 definitive secular variation model 
coefficients generated in step a. 

2.1 secular Variation Dat a Analysis (British Responsibility) 

The only data available for secular variation modeling are the 
observatory magnetic annual means, the first time derivative which 
provides information concerning the slow (greater than one year) rates 
of change of various components of the Earth , s main magnetic field 
at various geographic locations (roughly 200) around the world. 
Because of the sparsity and spatial nonuniformity of this data, it 
is possible to generate only a degree and order 8 spherical harmonic 
model of the secular variation. Furthermore, the predictive model 
is necessarily based on extrapolations of each magnetic component at 
each observatory site. Examples of observatory annual means from a 
few selected sites such as Honolulu, Huancayo, Pilar, and Rude Skov, 
for the X-north, Y-east, and Z-vertically down components of the 
Earth' s magnetic field, are given in figures (la) , (lb) , (lc) , (Id) , 
(le) and (If) through (4a), (4b), (4c), (4d) , (4e) and (4f) . The 
discontinuities in the field components at Honolulu are due to 
repositioning of the observatory at two separate instances. In several 
instances, the rate of change of one or more field components at an 
observatory has reversed direction over time intervals as short as 
two or three years. The sudden, unpredictable nature of the Earth's 
field is well illustrated by these observatories. The first-order 
time derivative of these data contains magnetic field contributions 
from the Earth's core as well as from the ionosphere and magnetosphere. 
It is difficult to remove the external field effects from these data 
because much of it is related to the solar cycle and many observatories 
do not have a sufficiently long history for a detailed analysis. 
Consequently, some external field effects are not removed from these 
data at the expense of a somewhat larger uncertainty in the secular 
variation model coefficients. 
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FIGURE la. NORTH X COMPONENT AT HONOLULU (HON). 
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FIGURE lb. EAST Y COMPONENT AT HONOLULU (HON). 
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FIGURE 1c. VERTICAL Z COMPONENT AT HONOLULU (HON) 
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FICURE Id. DECLINATION D COMPONENT AT HONOLULU (HON) 
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FIGURE le. INCLINATION I COMPONENT AT HONOLULU (HON). 




26 


FIGURE If. TOTAL INTENSITY F COMPONENT AT HONOLULU (HON) 
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FIGURE 2a. NORTH X COMPONENT AT HUANCAYO (HU A). 
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FIGURE 2b. EAST Y COMPONENT AT HUANCAYO (HUA) 
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FIGURE 2c. VERTICAL Z COMPONENT AT HUANCAYO (HUA) 
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FIGURE 2d. DECLINATION D COMPONENT AT HUANCAYO (HUA). 
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FIGURE 2c. INCLINATION I COMPONENT AT HUANCAYO (HUA). 
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FIGURE 3a. NORTH X COMPONENT AT PILAR (PIL). 
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FIGURE 3b. EAST Y COMPONENT AT PILAR (PIL). 
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FIGURE 3e. INCLINATION I COMPONENT AT PILAR (PIL) 
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FIGURE 3f. TOTAL INTENSITY F COMPONENT AT PILAR (P1L). 
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FIGURE 4a. NORTH X COMPONENT AT RUDE SKOV (RSV). 
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FIGURE 4b. EAST Y COMPONENT AT RUDE SKOV (RSV) 
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FIGURE 4c. VERTICAL Z COMPONENT AT RUDE SKOV (RSV). 
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FIGURE 4e. INCLINATION I COMPONENT AT RUDE SKOV (RSV). 
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FIGURE 4f. TOTAL INTENSITY F COMPONENT AT RUDE SKOV (RSV). 



Three secular variation models are generated by performing a weighted, 
least-square fit of the degree and order 8 spherical harmonic model 
to the first-order time derivative of the observatory annual means. 
These models were supplied by the British Geo log lea - urvey and are 
given in table 4 . 

2.2 Main ^ield Data A nalysis. (United States Respo..„xO_lity) 

The observatory annual magnetic means were not used in the main field 
modeling because those data contain, in addition to small extern-ij. 
field contributions, some rather large local and regional magnetic 
biases of crustal origin. A detailed survey at each observatory site 
would be necessary to remove these biases. Such surveys have rarely 
been performed due to the prohibitive cost, logistics, and 

international politics involved. 

The MAGSAT data consisted of 30,473 vector magnetic field values 
selected from 401 of the first 804 orbits. To minimize solar influences, 
the K/. magnetic index was required to be equal to, or less than, 
2(K/><2). These orbits were individually edited by an interactive 
graphics process to delete field aligned current effects and spurious 
data. Also, the following corrections for magnetospheric effects due 
to the ring current, magnetopause currents, and magnetotail currents 
were applied: 

B x (r,e,4>,T) = - < 7 °(x) sin 0 + {< 7 ,* (x) cos <|) + s, 1 (x) sin <J>} cos 0 (21a) 

B y (r,0,<}>,x) = sin <)>- s}(l) cos <]> (21b) 

B z (r,0,<|>,x) = < 7 ,°(x) cos 6 + { 4, 1 (x) cos <j> + s’ (x) sin <j)} sin 0 (21c) 


where the time-dependent coefficients are functions of the 
Storm Time (Dst) index: 

< 7 ,°(t) = 19.69 -0.63Dsr(x) 

< 7 j(x) = -0.38-0.06Dsr(x) 

5 i (x) = -2.90 + 0.17Ds:(x) 


Disturbance 


(22a) 

(22b) 

(22c) 


These corrections are derived from the external magnetic field 
potential : 
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TABLE 4. SECULAR VARIATION MODELS (units: nanotesles/year) 


1982.5 Epoch 1987.5 Epoch 1992.5 Epoch 


n m 


1 0 22.601 

1 10.491 

0 -14.454 

1 3.389 

2 5.043 

0 2.839 

1 -5.747 

2 -1.857 

3 -1.130 

0 0.000 

4 1 0.000 

4 2 -6.929 

4 3 0.000 

4 4 -6.192 

5 0 0.952 

5 1 -0.577 

5 2 -1.578 

5 3 -3.916 

5 4 0.000 

5 5 1.018 

6 0 0.962 

6 1 0.000 

6 2 1.678 

6 3 0.755 

6 4 0.000 

6 5 0.000 

6 6 1.437 

7 0 0.376 

7 1 -0.645 

7 2 0.457 

7 3 1.020 

7 4 1.583 

7 5 0.893 

7 6 0.387 

7 7 0.000 

8 0 0.651 

8 1 0.000 

8 2 -0.176 

8 3 0.000 

8 4 -0.842 

8 5 -0.222 

8 6 0.000 

8 7 -0.394 

8 8 -0.516 


0.000 18.745 0. 

-20.091 10.554 -17. 

0.000 -12.552 0. 

-14.476 3.290 -15. 

-20.686 0.656 -15. 

0.000 3.644 0. 

4.810 -6.918 4. 

3.013 0.000 1. 

-10.049 -4.759 -10. 

0.000 0.000 0. 

5.277 0.483 3. 

1.770 -7.391 2. 

4.323 0.510 3. 

0.771 -5.343 0, 

0.000 0.808 0, 

0.000 -0.357 -0 

-0.402 -1.730 0 

-0.501 -3.250 0 

0.000 0.000 1 

0.532 2.001 0 

0.000 1.296 0 

-1.107 0.000 0 

-0.821 1.799 -1 

-0.428 0.834 0 

-0.954 -0.667 -1 

0.349 0.000 0 

0.456 0.663 1 

0.000 0.670 0 

0.223 -0.483 0 

0.321 0.000 0 

0.375 1.010 0 

0.880 1.903 0 

0.000 0.597 0 

0.467 0.000 0 

0.825 0.543 0 

0.000 0.237 0 

0.459 -0.692 0 

-0.332 -0.367 -0 

0.330 0.000 0 

-0.268 -1.292 0 

0.428 0.000 0 

-0.967 0.423 -0 

-0.994 0.000 -0 

0.693 -0.562 0. 


000 16.013 0.000 

836 9.259 -13.759 

000 -11.703 0.000 

523 3.715 -12.790 

631 1.767 -14.865 

000 2.115 0.000 

285 -7.596 3.082 

498 0.000 0.844 

351 -5.815 -11.342 

000 -0.770 0.000 

375 0.968 3.281 

994 -7.414 3.680 

812 0.775 2.799 

000 -6.361 0.000 

000 1.662 0.000 
.593 0.000 0.000 
.000 0.000 -2.096 

.000 -2.699 1.226 

.712 0.000 1.193 

.464 3.001 0.650 

.000 0.751 0.000 

.000 0.000 -0.583 

.277 1.451 -0.644 

.000 0.000 0.000 

.664 0.000 -2.266 

.000 0.000 0.000 

.381 0.000 0.000 

.000 0.505 0.000 

.978 0.000 0.599 

.000 -0.869 0.793 

.390 1.457 0.000 

.000 2.650 0.000 

.000 - 1.020 0.000 

.000 0.000 0.417 

.000 0.000 0.000 

.000 0.000 0.000 

.399 -1.089 0.427 

.195 0.000 -0.809 

.586 0.000 0.507 

.519 -2.114 0.349 

.528 0.000 0.499 

.742 0.978 0.000 

.713 0.000 -0.684 

0000 0.000 0.000 
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' V cu" « f r V 




{q*(x) cos m <J) + 5 *(t) sin mty }PT(cos 0 ) 


(23) 


when N txt = 1, via the relations: 

i ay 

B * = “Be “ r 


(24a) 


B v 


= +B* - 


1 


dV 


rsin 0 3 <J> 


(24) 


dV 

B z - -B. - 37 


(24c) 


A further correction takes into account the magnetic fields induced 
in the Earth by the external fields which, because of their time 
dependence and the generally low but finite conductivity of the crust 
and mantle, induce electric currents in the crust and mantle, which 
in turn generate secondary magnetic fields. These secondary fields 

are of internal origin and primarily affect the g° coefficients of 
the internal magnetic potential: 


N <* "fa Y + ' 

V te (r,e,fc X) = a I 21 - 

n m 1 m * 0 \^ r J 


{ g?(*) COS mfy + h” (t) sin m <|> }PT (cos 0 ) 


; 25) 


Taking derivatives as before with n -1 and m- 0 , the magnetic field 
corrections due to induction effects are: 


B x (r, 0 , 4 >,T) = - 

B r (r,0,<M) = 0 


f-T 


£,°iC0 sin0 


fa Y 

B z (r, 0 ,<M) = - 2 1 -J«„(x)cos0 

where the induced part of the g? coefficient is given as: 


(26a) 

(26b) 

(26c) 


g\lO = 0.27tf(t) 
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The external and induced magnetic field corrections given above were 
subtracted from the MAGSAT observations. These corrections are based 
on previous analyses of MAGSAT data by Langel and Estes (1985) and 
by Quinn, Kerridge, and Barraclough (1986) . 


No attempts were made to remove magnetic influences due to ionospheric 
currents such as those generated by solar quiet (SQ ) currents, auroral 
electro jet currents, and equatorial electrojet currents, which are 
located below the MAGSAT orbit altitudes. These influences, thong., 
generated external to the Earth's surface, are nevertheless part or 
because their sources are internal to the point of 
observation. Consequently, separating core-generated fields from 
crustal and ionospherically generated fields measured by satellite 
magnetometers is difficult and is still a research matter . Fortunately, 
fields generated in the Earth's crust and ionosphere are significantly 
attenuated at satellite altitudes. Therefore, errors in the main 
field model coefficients due to contamination of the satellite data 
by these two sources are comparatively small. 


MAGSAT was 
(NASA) /U . S . 
NASA in the 


a joint National Aeronautics and Space Administration 
Geological Survey mission. These data were supplied by 
form of investigator B tapes. 


The DE-2 satellite data set consisted of 5,100 data points gleaned 
from the low-altitude end of a comparatively eccentric orbit. is 
data set contained only scalar total intensity measurements of the 
Forth's maanetic field. However, this data set exhibited 
substantially higher rms error relative to the DGRf-80 model than the 
MAGSAT data, consequently, for this and other reasons, the DE-2 data 
were not used in the final main field model determination. This data 
set originated with M. Sugiura of Japan (formerly of NASA) . It was 
edited by J.R. Ridgeway of Science Applications Research Corporation 
and it was subsequently supplied to NAVOCEANO by Dr. Robert Langel 

Of NASA. 


The Project MAGNET aeromagnetic data consisted of 338 high level 
flights (^ 15,000 feet) of vector component measurements. These data 
are routinely processed by NAVOCEANO at a 2-second sample rate and 
sent to the National Geophysical Data Center (NGDC) in Boulder 
Colorado. A weak low-pass filter with a cut off wavelength of 
approximately 7 km is routinely applied to this high-level data. The 
cut-off wavelength will vary slightly, depending on the average speed 
of the aircraft, which depends on prevailing wind conditions at the 
time of flight. Typical flights last 10 to 12 hours and are flown 
at an av age -speed of 440 km/hr. They are generally flown at r.ight 
in orde to minimize solar-driven external field effects which 

contribute to the Daily Variation ( DV ) of the Earth's field. Project 
MAGNET flights are of long range in remote ocean areas, which precludes 
the monitoring of DV . Therefore, no explicit DV corrections are 
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made to the data. Also, the aircraft's vector magnetometer is 
calibrated at the NASA Coil Room Facility at the Goddard Space Flight 
Center in Maryland at least once a year. 

The magnetic field observations returned from each Project MAGNET 
flight are routinely reduced in accordance with the following 
procedures : 

a. Rotate vector measurements from magnetometer coordinates to 
instantaneous aircraft coordinates. This rotation involves only sma.l 
misalignments relative to an imaginary coordinate system rigidly 
attached to the aircraft. 

b. Compensate in aircraft coordinates for the perturbing magnetic 
effects associated with the presence of the aircraft by removing a 
field phenomenologically modeled as: 


B c (t) 


B 


Perm 


+ a B m (t) 


+ P 


d B„(t) 

dx 


( 29 ) 


where the first term represents the permanent magnetic field generated 
by the remnant magnetization of the aircraft's metal parts, the second 
term represents the field induced in the aircraft's metal structure 

by the presence of the ambient field B m (t), and the third term represents 
magnetic fields generated by eddy currents created on the aircraft's 
metal surfaces by the aircraft's motion through the Earth's spatially 
varying field. Here, B m (t) is the magnetic field measured by the 
magnetometer after it has been rotated into aircraft coordinates as 
indicated in step a. The compensation model contains 21 coefficients, 

3 in the vector B Ptrm , 9 in the 3x3 matrix a, and 9 in the 3x3 matrix 

P- 


c. Rotate the compensated field from instantaneous aircraft 
coordinates to geodetic coordinates, taking into account the 
misalignment of the inertial attitude device relative to the 
instantaneous aircraft coordinates. 

d. Visually edit the data in each flight via interactive graphics 
techniques . 

The compensation coefficients are determined by performing calibration 
flights at an altitude of 1,500 feet above a designated magnetic 
observatory. “These flights consist of a set of yaw, pitch, and roll 
maneuvers performed along the four cardinal headings (north, south, 
east, and west) . The coefficients are then determined by a 
least-squares procedure that minimizes the squared difference between 
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fnnwArd continued and rotated into instantaneous 

attitude devices on the 

aircraft) and the field measured by the aircraft's magnetometer. 

Using this minimization technique, the^^compensation coefficients 
are determined simultaneous y t a anc ^ 3 f or the inertial 

magnetometer misalignment A s tep c The overall process 

attitude device misalignment mentioned in St. iterations to converge, 
is, therefore, nonUnear requiring the relati-a 

Note, however, that and the inertial system 

misalignment between the magneto angles are arbitrarily 

-fto L h ?o r % f o 0 th'ae\n th p?a=tTc?, e riythe three inertial system bias 
angles (Euler angles) are computed. 

. . _ nr i editing, the aeromagnetic data were decimated 
After compensation and edl r q ' every hundredth point was 

to a 200-second sample i“erval 'Ve t ic observations . The resulting 
selected) , yielding 54, 656 vector inverted from geodetic coordinates 

to°s j p e he” E Ltrd a inltes using the coordinate transformations of the 

previous section. 

For modeling, it is desl |^^ t0 Consequently^th^MAGSAT data set, 
or backward to a common p * roor dinates, was pushed forward to 

which was originally s^s^Sfar variation model. The portion of the 
1985.0 via the i982 . 5 sec ui .ar vet . l 1985.0 was pushed forward 

Project 1982 5 secular variation model, while the Portion 

to 1985.0 via Jmft data set collected after 1985.0 was pushed 

of the Project ^?!^L1 q87 I secular variation model. Subsequently, 
backward to 1985 . 0 by the 198 . , , g aenerated by performing 

a revised 1985.0 epoch mam f: ie d Spherical 

a weighted least-squares fit o and Project MAGNET data sets. 

harmonic model to the combined ^GSAT^d Project 19g? _ 5 secular 

The resulting main field m0 ® ' WC-85 (revised) . These coefficients 

variation model, is referred c in fie xd model was produced by 

are listed in table 5. The 1 . 99 °;"” a '-"“xcal harmonic coefficients 

forward fnNTm'e us'in^ the i M7?S .«»Ur ' variation =oem^en« The 

r^SciJr.-^orSJi’iS coefficients of 

which are listed in table 3 of section 1.3. 

2.3 Derails of Ma i n FI M '1 Hode li a a 

Cain^t 6 ^!? 5 ^gglf^^Th^tjlTj^clfive^as't^minimiz^'the r chi-square 

function 

(29) 


(x 2 ) 


X 2 = X? + Xe + xJ + Xf 
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TABLE 5. WC-85 (REVISED) SCHMIDT NORMALIZED GAUSS COEFFICIENTS 


n 

m 

s ” 

( nT ) 

K 

( nT ) 

c 

( nT / yr ) 

K 

( nT / yr ) 

i 

0 

- 29874.2 

.0 

18.7 

. 0 

i 

1 

- 1904.5 

5496 . 4 

10 . 6 

- 17 . 8 

2 

0 

- 2071 . 6 

. 0 

- 12.6 

. 0 

2 

1 

3045.7 

- 2200 . 6 

3.3 

- 15.5 

2 

2 

1688.7 

- 306.1 

.7 

-15 . 6 

3 

0 

1294.7 

. 0 

3 . 6 

.0 

3 

1 

- 2210 . 1 

- 306.4 

-6 . 9 

4 . 3 

3 

2 

1246.8 

284.2 

.0 

1 . 5 

3 

3 

832 . 4 

- 300.7 

-4 . 8 

- 10 . 4 

4 

0 

933.5 

.0 

. 0 

.0 

4 

1 

782.5 

232.5 

.5 

3.4 

4 

2 

360.5 

- 247 . 6 

- 7.4 

3.0 

4 

3 

- 424.2 

72.2 

.5 

3.8 

4 

4 

166.0 

- 296.5 

- 5.3 

.0 

5 

0 

- 212 . 3 

.0 

. 8 

.0 

5 

1 

354 . 0 

43.7 

-.4 

-.6 

5 

2 

255.2 

148.7 

- 1.7 

.0 

5 

3 

- 94 . 6 

- 154 . 6 

-3 . 3 

.0 

5 

4 

- 162.3 

- 76.2 

.0 

1 .7 

5 

5 

- 47.2 

95.0 

2.0 

.5 

6 

0 

52.5 

.0 

1 . 3 

.0 

6 

1 

63.7 

- 14.7 

. 0 

.0 

6 

2 

51.0 

88.6 

1.8 

-1 . 3 

6 

3 

- 185.4 

70.0 

.8 

.0 

6 

4 

3.8 

- 47.8 

-.7 

- 1.7 

6 

5 

15.4 

- 1.4 

.0 

.0 

6 

6 

- 99.3 

17.7 

.7 

1 . 4 

7 

0 

72.8 

.0 

.7 

.0 

7 

1 

- 59.7 

- 83.5 

-.5 

1 . 0 

7 

2 

1.3 

- 26.7 

.0 

.0 

7 

3 

25.1 

-1 . 9 

1.0 

.4 

7 

4 

- 4.8 

19.9 

1 . 9 

. 0 

7 

5 

4.9 

17.9 

. 6 

.0 

7 

6 

10 . 1 

- 21.5 

.0 

. 0 

7 

7 

-.8 

- 6.8 

.5 

. 0 
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TABLE 5. WC-85 (REVISED) SCHMIDT NORMALIZED GAUSS COEFFICIENTS (con.) 
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TABLE 5. WC-85 (REVISED) SCHMIDT NORMALIZED GAUSS COEFFICIENTS (con.) 


n 

m 

m 

gn 

( nT ) 

i m 

K 

( nT ) 

. m 

gn 

( nT / yr ) 

hn 

( nT / yr ) 

11 

0 

1 . 3 

.0 

.0 

.0 

11 

1 

-1 . 4 

.0 

.0 

.0 

ii 

2 

- 2.5 

1.0 

.0 

. 0 

ii 

3 

3.2 

- 1.6 

.0 

. 0 

ii 

4 

.2 

- 2.2 

.0 

.0 

ii 

5 

-1 . 1 

1 . 1 

.0 

. 0 

ii 

6 

.3 

-.7 

.0 

.0 

ii 

7 

-.3 

- 1.7 

.0 

.0 

ii 

8 

. 9 

- 1.5 

.0 

.0 

ii 

9 

-1 . 1 

- 1.3 

. 0 

.0 

ii 

10 

2 . 4 

-1 . 1 

.0 

.0 

ii 

11 

3.0 

.6 

.0 

. 0 

12 

0 

- 1.3 

.0 

.0 

. 0 

12 

1 

.1 

.7 

.0 

.0 

12 

2 

.5 

.7 

.0 

.0 

12 

3 

.7 

1.3 

.0 

. 0 

12 

4 

.4 

- 1.5 

.0 

. 0 

12 

5 

-.2 

.3 

.0 

. 0 

12 

6 

- 1.1 

.2 

.0 

. 0 

12 

7 

.9 

- 1.1 

.0 

.0 

12 

8 

-.6 

1.2 

.0 

.0 

12 

9 

.8 

-.2 

.0 

. 0 

12 

10 

.2 

- 1.3 

.0 

. 0 

12 

11 

.4 

.6 

.0 

.0 

12 

12 

.2 

. 6 

.0 

. 0 


S3 










with respect to the '168 internal Gauss coefficients of a degree and 
order 12 spherical harmonic model, where: 


i, 


'"i ru 

II 

2>„(B ,,-bJ 
1 = 1 

II 

C4 

-o 

1 

CQ 

& 

* 

s ^*4 « 

II 

Ivv„(B„-&*) 2 

i = 1 

II 

V 

i, w fl (B Fl -b F ,y 

t = 1 


(30a) 


(30b) 


(30c) 


( 30d) 


where the upper case B' s refer to the model values of their respective 
magnetic components, while the lower case b's refer to the observed 
(measured) values of their respective magnetic components. Th 
subscript i refers to a particular data point, the total number I of 
which may differ for each magnetic component. Each data point is 
weight e d ^b y a weight factor, w, which depends on several factors: 


a. Data type W m 


MAG SAT = 1 


Project MAGNET — 1/4 


Project MAGNET observatory airswing calibrations yield rms errors on 
the order of 35 nT, while MAGSAT rms differences from degree 12 
spherical harmonic models yield rms values °n the order of 9 nT 
Consequently, the relative weight of the two data sets is 

b ~l~l T his factor characterizes the relative quality of the two 
data sets. 


b The relative number of data points per equal area (5 x5 at 
the equator) cell; each cell was given equal weight. Therefore, dat 
points corresponding to cells with more than the average number 

points per cell, N, received less weight and vice versa. 


c The relative rms error of data in a particular flight or orbit 
relative to the rms error, O, for all data of the corresponding data 
type (MAGSAT or Project MAGNET) . 
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d. The relative rms error of data of a specified type in an equal 
area cell relative to all data of that type, o. 

e. The age of the data relative to the model epoch 1985.0. Thus^ 
data collected five years away from this epoch get a weight of 
approximately 1/3, while data collected at the model epoch get a 
weight of 1 . 


f . 

Distance of geomagnetic latitude, ©, f 

equator 

• 





ie„ 

|< 20° | 

= 1 
= 0 

for 

for 

k = l,2,3(r,e,4>)] 
k=4(F) J 

| ; n = 1 

ie„ 

j> 20° | 

= 0 
= 1 

for 

for 

k = 1,2, 3(r,0,<t>)l 
k = 4(F) 

| ; n= 2 


Ola) 


(31b) 


This weighting scheme then takes the following mathematical form. 


WiJUm. = 


( jsj \( rr 


km 


hm Ukm 

®kmi J v ®kmj 


(?T 


(32 


where the indices correspond to the following: 
ith - data point 

jth - equal area cell (1654 total) 
kth - magnetic component (r,0,<j>,F) 

Ith - aircraft flight or satellite orbit 
mth - data type (MAGSAT, Project MAGNET) 
nth - geomagnetic latitude band (n = 1 or n- 2) 

The decay constant X was arbitrarily chosen to be 5 years, while 
At, = X, — T epoch 


where x, is the time of observation in years and T epoch xs 1985.0. 

Table 6 gives the overall rms errors of a particular magnetic component 
for each of the three separate data sets relative to the DGRF/IGRF 
series of SIMMs'. Table 7 lists the number of data points associated 
with each mangetic component for each of the three data sets. Table 
8 lists the average number of data points per 5 x5 equal area cell 
for each magnetic component for each data set. Rms statistics relative 
to the DGRF/IGRF series of models for the Project MAGNET data set are 
further broken down by Project ID and flight number in table 9. Due 
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TABLE 6. RMS ERRORS RELATIVE TO IGRF/DGRF MODELS 



rms (nT) 

rms (nT) 

rms (nT) 

rms (nT) 

MAG SAT 

14.7 

12.0 

13.2 

12.7 

Project MAGNET 

101.7 

107 . 9 

105.5 

96 . 4 

DE-2 

— 

1 

1 

1 



1 

1 

1 

122 . 4 


TABLE 7 . NUMBER OF RECORDS 



N. 


N, 

n f 

MAGSAT 

30473 

30473 

30473 

30473 

Project MAGNET 

54656 

54656 

54656 

54656 

DE-2 

1 

1 

1 

1 

1 

1 

I 

1 

5100 


TABLE 8 . AVERAGE NUMBER OF RECORDS PER CELL 



n m 

Ny 

N, 

N, 

MAGSAT 

18.4 

18.4 

18.4 

18 . 4 

Project MAGNET 

33.0 

33.0 

33.0 

33.0 

DE-2 

1 

1 

1 

1 

1 

1 

— 

3.1 
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TABLES 9. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO ICRF/DGRF 
MODELS (RMS units: nT) 


PROJECT 

FLIGHT 

A32-153 

6005 

A32-153 

6006 

A32-153 

6007 

A32-153 

6008 

A32-153 

6009 

A32-153 

6010 

A32-153 

6011 

A32-153 

6012 

A32-153 

6013 

A32-153 

6014 

A32-153 

6015 

C32-252 

1074 

C32-252 

2007 

C32-252 

2009 

C32-252 

2011 

C32-252 

2012 

C32-252 

2013 

C32-252 

2014 

C32-253 

1077 

C32-253 

1080 

C32-253 

1081 

C32-253 

1084 

C32-253 

1085 

C32-253 

1086 

C32-253 

1087 

C32-253 

1088 

C32-253 

1091 

C32-254 

4027 

C23-245 

4028 

C32-254 

4049 

C32-254 

5022 

C32-254 

5023 

C32-254 

5026 

C32-254 

5037 

C32-254 

5030 

C32-254 

5031 

C32-351 

1097 

C32-351 

3074 

C32-351 

30"’ " 

C32-351 

405 . 

C32-351 

4032 

C32-351 

4033 

C32-351 

4034 

C32-351 

4035 

C32-351 

4036 

C32-351 

5033 



JULIAN 

RMS 

YEAR 

EM 

X 

1981 

129 

116.4 

1981 

133 

141.4 

1981 

139 

223.7 

1981 

143 

88.5 

1981 

146 

138.3 

1981 

149 

132.5 

1981 

153 

200 . 1 

1981 

155 

107.8 

1981 

157 

283.0 

1981 

159 

146 . 4 

1981 

162 

111.3 

1982 

62 

60 . 1 

1982 

64 

68 . 9 

1982 

67 

79.4 

1982 

74 

84 . 0 

1982 

77 

95.1 

1982 

89 

63.4 

1982 

96 

80.5 

1982 

135 

85.4 

1982 

149 

86.0 

1982 

161 

133.6 

1982 

171 

63.0 

1982 

175 

58.4 

1982 

177 

164.8 

1982 

180 

98.2 

1982 

182 

122.5 

1982 

186 

112.4 

1982 

231 

86.4 

1982 

233 

155.3 

1982 

236 

120 . 9 

1982 

243 

100.2 

1982 

246 

87 . 9 

1982 

251 

88.4 

1982 

260 

81.1 

1982 

274 

81 . 1 

1982 

277 

131.6 

1982 

304 

150 . 6 

1982 

305 

109.4 

1982 

307 

105.5 

1982 

309 

103 . 6 

1982 

312 

86.5 

1982 

314 

94 . 9 

1982 

320 

73 . 0 

1982 

322 

59 . 1 

1982 

324 

71 . 6 

1982 

327 

208.2 


RMS 

RMS 

RMS 

1 

z. 

E 

92 . 3 

159.0 

156 . 6 

98.4 

106.1 

104 . 6 

242 . 1 

191.5 

200 . 9 

121 .4 

138 . 6 

137 . 5 

86.6 

105 . 3 

113.9 

135.3 

151 . 6 

151.5 

115.8 

17" . 9 

176 . 0 

233.3 

107 . 0 

106 . 8 

128 . 6 

181.6 

175.2 

166.2 

110.6 

111 . 3 

101 .2 

135.2 

135 . 8 

51.9 

49.2 

73.2 

96.6 

87 .4 

82 . 1 

76 . 6 

85.1 

74 . 0 

108.4 

103 . 9 

103.4 

74 . 4 

65.7 

94 .5 

60.0 

64 . 8 

61 . 8 

70.8 

79.8 

76 . 8 

68.7 

93.7 

93 . 8 

121 . 8 

58 . 9 

66 . 3 

48.8 

48 . 8 

126.9 

54.8 

99 . 1 

91 . 1 

61 . 6 

78 . 1 

73.5 

103.3 

162.0 

185 . 9 

87 . 1 

119.7 

122.7 

134.2 

160.7 

154 . 8 

66.2 

98 . 3 

100 . 1 

91.4 

64 . 1 

70.2 

81.4 

81 . 1 

16. 4 

87.9 

89.2 

78 . 1 

71.1 

50.6 

60 . 9 

74.2 

67.2 

60 . 8 

63.5 

71 . 6 

84.4 

79.5 

60.3 

83.0 

80.3 

8 3 . 0 

84 . 0 

46.9 

90.3 

81 . 1 

91.3 

128.3 

147 . 6 

126.3 

63.6 

101.4 

137.6 

89 . 6 

99 . 9 

57.3 

79.9 

96 . 1 

128.3 

88 . 9 

91.8 

101 . 6 

97 . 9 

100.0 

92.1 

115.1 

93 . 1 

144 . 1 

64 . 6 

66.8 

85.6 

105.2 

100 . 3 

191.1 

118.5 

124 . 3 
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TABLES 9. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO IGRF/DGRF 
MODELS (RMS units: nT)(con.) 


PROJECT 

FLIGHT 

YEAR 

JULIAN 

DAY 

C32-351 

4037 

1982 

332 

C32-351 

4038 

1982 

334 

C32-351 

4039 

1982 

337 

C32-351 

4040 

1982 

339 

C32-351 

4041 

1982 

341 

C32-351 

4042 

1982 

345 

C32-351 

4044 

1982 

350 

C32-351 

1098 

1982 

351 

C32-352 

1089 

1983 

27 

C32-352 

1090 

1983 

31 

C32-352 

1099 

1983 

34 

C32-352 

5045 

1983 

75 

C32-352 

5046 

1983 

78 

C32-352 

5047 

1983 

81 

C32-352 

5048 

1983 

84 

C32-352 

4045 

1983 

87 

C32-352 

4046 

1983 

91 

C32-352 

4048 

1983 

96 

C32-352 

1100 

1983 

99 

C32-353 

3076 

1983 

228 

C32-353 

3077 

1983 

130 

C32-353 

3080 

1983 

161 

C32-353 

3081 

1983 

166 

C32-353 

3082 

1983 

169 

C32-353 

3083 

1983 

171 

C32-353 

3084 

1983 

175 

C32-353 

3085 

1983 

179 

C32-353 

3086 

1983 

183 

C32-353 

3087 

1983 

186 

C32-451 

3088 

1983 

325 

C32-451 

3089 

1983 

327 

C32-451 

3090 

1983 

330 

C32-451 

5050 

1983 

144 

C32-451 

5051 

1983 

336 

C32-451 

4051 

1983 

340 

C32-451 

4052 

1983 

343 

C32-451 

4053 

1983 

346 

C32-451 

4054 

1983 

349 

C32-352 

4046 

1983 

91 

C32-451 

1103 

1983 

354 

C32-452 

1108 

1984 

35 

C32-452 

1106 

1984 

41 

C32-452 

-504 9 

1984 

45 

C32-452 

3091 

1984 

51 

C32-452 

5052 

1984 

64 

C32-452 

4050 

1984 

67 

C32-452 

4055 

1984 

73 

C32-452 

4057 

1984 

76 

C32-452 

4058 

1984 

78 

C32-452 

4059 

1984 

84 

C32-452 

4060 

1984 

88 


RMS 

RMS 

RMS 

RMS 

X 

1 

Z 

E 

64.3 

75.7 

69.8 

69.5 

83.3 

87 . 6 

89.7 

76 . 5 

72 . 9 

91 . 8 

90 . 8 

84 . 5 

52 . 6 

67 . 9 

90 . 9 

90.2 

146.6 

104 . 0 

81.2 

135 . 6 

83.8 

78 . 9 

52 . 6 

82 . 9 

63.0 

83.2 

86 . 0 

73.5 

129.4 

67 . 8 

145 . 8 

160 . 9 

64.4 

65.0 

81 . 3 

78 . 5 

50.1 

80.2 

37 . 3 

44 . 1 

75 . 6 

75.2 

101 . 3 

87.3 

98.9 

101.5 

79 . 9 

98.0 

89.6 

80.5 

81.6 

80.3 

109.5 

72.7 

83.2 

75 . 3 

212.6 

106.5 

145.5 

130.5 

255.4 

236.2 

116.9 

138.8 

39.4 

76 . 6 

56.8 

45.1 

59.2 

93.9 

58.7 

69.6 

31.8 

35.3 

71 . 1 

64 . 3 

86.3 

76.8 

111.5 

114 . 6 

64 . 1 

191.7 

94 . 8 

76.7 

91.0 

71 . 4 

67 . 8 

94 . 8 

124.0 

61.9 

115.2 

136 . 6 

96.1 

51.1 

98.2 

93.4 

59.7 

121.5 

76.2 

67 . 0 

73.2 

110.1 

96.5 

93.4 

91.6 

100.1 

103 . 9 

114.5 

82.7 

64.7 

87.7 

105 . 9 

44.8 

100 . 9 

64.8 

58 . 9 

73.4 

73.7 

78.4 

83 . 8 

92.5 

66.1 

60 . 9 

96 . 1 

107.6 

169.3 

112.0 

125.5 

93.4 

81 . 6 

138.8 

134 . 3 

48.9 

49 . 9 

56.8 

57 . 6 

174.0 

133.1 

103.2 

114 . 6 

114 . 6 

67.4 

95.5 

105.8 

109.1 

130.4 

109.4 

107.4 

77.8 

74.0 

94.5 

101 . 6 

39.4 

7 6.6 

56.8 

45 . 1 

29.7 

35 . 8 

52 . 8 

47.5 

103.7 

86.2 

110 . 9 

54 . 4 

115.3 

113.8 

124 . 9 

122 . 6 

105.9 

58.8 

156.2 

63.5 

53.2 

265.5 

242 . 9 

56.3 

120.5 

194.8 

111 . 9 

97.0 

168.6 

175.2 

80.7 

64.3 

146.1 

91 . 9 

174.5 

101.3 

61.8 

58 . 1 

123.0 

44.3 

86.2 

68 . 6 

127.7 

55.3 

51.3 

139.6 

158 . 6 

45.2 

28.8 

138.5 

98.0 

38 . 3 


58 


TABLE 9. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO ICRF/DCRF 
MODELS (RMS units: nT) (con.) 





JULIAN 

RMS 

RMS 

RMS 

*MS 

PROJECT 

FLIGHT 

YEAR 

DAY 

X 

1 

z. 

p 

C32-452 

4061 

1984 

91 

76.2 

62.7 

93. 1 

90 . 9 

C32-452 

1107 

1984 

93 

105,8 

79.1 

139.6 

112.5 

C32-453 

1109 

1984 

137 

97.5 

86.3 

116.8 

115.4 

C32-453 

3092 

1984 

140 

163.7 

83.1 

96.5 

83 . 5 

C32-453 

3094 

1984 

157 

80 . 1 

99 . 9 

119.2 

86 . 9 

C32-453 

3095 

1984 

161 

70.8 

64 . 8 

63.8 

79 . 6 

C32-453 

3096 

1984 

164 

71.5 

153.6 

121 . 6 

68 .7 

C32-453 

3097 

1984 

171 

100.8 

71.3 

88.8 

69 . 9 

C32-453 

3098 

1984 

168 

138.9 

149.5 

115.5 

81 . 9 

C32-453 

3099 

1984 

171 

114.3 

79.5 

98 . 6 

107.2 

C32-453 

3101 

1984 

176 

61.4 

72.5 

75 . 3 

76.2 

C32-453 

3100 

1984 

180 

98.4 

83.6 

113 . 5 

118.0 

C32-453 

3102 

1984 

183 

78.8 

69.8 

101 .2 

101 . 8 

C32-453 

3103 

1984 

186 

96.7 

78.4 

125 . 3 

113.1 

C32-453 

3104 

1984 

190 

68.7 

84 . 4 

79.0 

88.2 

C32-453 

1110 

1984 

192 

83.5 

88 . 6 

127.5 

122 . 6 

C32-454 

1048 

1984 

215 

135.4 

111.5 

165.6 

165.2 

C32-454 

1075 

1984 

219 

62.7 

61.5 

96.2 

96 . 6 

C32-454 

1070 

1984 

221 

59.2 

74.2 

96 . 3 

95 . 8 

C32-454 

1121 

1984 

224 

75.6 

86.2 

100 . 3 

84.0 

C32-454 

1122 

1984 

227 

69.4 

71.7 

88.2 

70.5 

C32-454 

1123 

1984 

230 

54.4 

52.2 

77 . 8 

70 . 1 

C32-454 

1124 

1984 

237 

59.7 

62.5 

100 . 3 

106 . 8 

C32-454 

1117 

1984 

242 

101.4 

79.2 

126 . 6 

129 . 0 

C32-454 

1118 

1984 

246 

44.2 

78.5 

66.4 

57.9 

C32-454 

1119 

1984 

249 

52 . 6 

55.5 

71 . 1 

61 . 1 

C32-454 

1120 

1984 

252 

57.9 

41.8 

83 . 9 

66 . 1 

C32-454 

1111 

1984 

255 

74.4 

63.4 

101 . 4 

100 . 0 

C32-454 

1112 

1984 

258 

93.6 

91.6 

113.1 

111 . 5 

C32-454 

1114 

1984 

262 

66.3 

60.8 

42 . 3 

55.4 

C32-454 

1115 

1984 

272 

34.8 

37.4 

57.2 

52.3 

C32-551 

1126 

1984 

290 

64.6 

69.8 

81 . 9 

70 . 9 

C32-551 

2015 

1984 

292 

117.5 

134.4 

44.2 

141.8 

C32-551 

2016 

1984 

295 

54.2 

101 . 1 

97 . 9 

94.2 

C32-551 

7005 

1984 

298 

174.7 

91.3 

114 . 3 

112.2 

C32-551 

7006 

1984 

304 

80.4 

97.5 

75.4 

76.0 

C32-551 

7007 

1984 

307 

218.4 

186.7 

96.3 

96.8 

C32-551 

4065 

1984 

310 

206.9 

218.9 

121 . 1 

85.7 

C32-551 

4066 

1984 

313 

101.5 

103.3 

81 . 6 

79.7 

C32-551 

4067 

1984 

319 

70.7 

137.9 

87 . 8 

59 . 1 

C32-551 

4068 

1984 

321 

198.6 

160.0 

107.2 

86 . 6 

C32-551 

4069 

1984 

325 

119.4 

170.7 

97.7 

62 . 9 

C32-551 

4070 

1984 

329 

76.6 

232.9 

79.9 

63 . 3 

C32-551 

4071 

1984 

333 

121.9 

124 . 8 

74 . 6 

69 . 4 

C32-551 

4072 

1984 

337 

91.2 

109.8 

104 . 4 

93.2 

C32-551 

3105 

1984 

342 

84.6 

52.4 

150.7 

97 . 3 

C32-551 

1127 

1984 

345 

45.9 

48.7 

95.5 

87 . 6 

C32-552 

1125 

1985 

16 

81.6 

61.5 

77 . 1 

54 . 4 

C32-552 

1131 

1985 

19 

87.2 

105.9 

85.2 

54 . 0 

C32-552 

1132 

1985 

27 

115.7 

106.3 

72.2 

62 . 7 

C32-552 

1133 

1985 

30 

88.6 

260.8 

138.9 

124 . 5 
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TABLE 9. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO iCRF/DGRF 
MODELS (RMS units: nT) (con.) 


PROJECT 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-552 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-553 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-554 

C32-651 

C32-651 

C32-651 

C32-651 

C32-651 

C32-651 

C32-651 


FLIGHT 

2017 

2029 

2030 

2031 

2032 

2033 

2034 
2022 
2023A 

2023 

7009 
7008 

7010 

7011 
2036 

1129 

1130 
1134 

3106 

3108 

3109 

3110 

3107 

3111 

3112 

3113 

3114 

3116 

3117 

3118 

3119 
3121 
3122 

1135 

1137 

1136 

1138 

1139 

1140 
6016 
1144 

1141 

1142 
6017 

1143 

4074 

4075 

3125 

3126 
5053 
4077 


YEAR 

1985 

1985 

1985 

1985 

1985 

1985 

1985 

1985 

1985 

1985 

1984 

1985 

1984 

1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985' 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 
1985 


JULIAN 

DAY 

33 

35 

37 

40 

43 

46 

49 

53 

57 

61 

62 

66 

69 

71 

75 

80 

84 

114 

116 

120 

123 

126 

129 

136 

140 

144 

148 

154 

157 

160 

169 

171 

174 

196 

210 

202 

213 

217 

220 

230 

244 

249 

250 
252 
281 
285 
290 
301 
305 
314 
317 


RMS 

X 

58.0 

46.1 

68 . 9 

85.4 

62.6 

54 . 4 

35.3 

33.6 

40.0 

75.3 
114.5 

73.0 

101.7 

51.4 

78 . 6 

209.1 

63.3 

135.9 

54.1 

57 . 3 

62.7 

59.2 

50.1 

101.9 

60.6 

93.0 

61.0 

81.2 

132.3 

129.8 

223.0 

232.4 

134.8 

105.3 

62.2 

60.8 

109.1 

68.0 

155.7 

151.3 

147.6 

96.4 

155.0 

182.3 

65.6 

32.3 

38.0 

72.0 

61.4 

77 . 3 

97.0 


RMS 

1 

136.3 

182 . 4 

198.5 

127.8 

137.2 

113.0 

110 . 4 

116 . 8 

125.1 

127.6 

99.4 

103.3 

104 . 1 

123.0 

158.4 

201 . 6 

59.0 

130.2 

112.2 

112.8 

112.7 

99.9 

131.1 

63.2 

79.0 

67.1 

158.4 

75.5 

82.7 

89.1 

145.9 

111.7 

193.1 

123.5 

58.0 

89.1 

115.5 

78.5 

109.2 

158.4 

116.2 

83.0 

128.9 

94.8 

101.2 

148.1 

225.2 

236.2 

110.8 

178.2 

202.1 


RMS 

Z. 

49.0 

41.7 

59 . 2 

46.0 

104.0 

73 . 8 
48 . 6 

59.5 
64 . 1 

89.7 

107.0 

75 . 8 

93.0 

65.5 

83.4 

127.2 

71.4 

96 . 9 

88.0 

74 . 9 

55.0 

59.2 

52 . 2 

157.6 

145 . 9 

82.7 

80 . 6 

120.8 

120 . 9 

145.6 

114.4 

127 . 1 

130.4 

137.1 

69.9 

85.9 

111.9 

85.3 

126.9 

165 . 6 

186.4 

91.3 

151.0 

136.2 

88.1 

83.1 

85.0 

95.8 

80.4 

105.1 

114.7 


RMS 

F 

65 . 1 

39.0 

78.5 

105.1 

111.7 

71.2 

40 . 4 

46.1 

51 . 9 

85.4 

132.0 

73.0 

102 . 1 

68.5 
83 . 6 

51.3 

70 . 1 

75 . 0 

92 . 6 

72.2 

56.2 

56.4 

55.8 

72 . 9 

61.2 

89.7 

61.3 

61.2 

75 . 1 

106.9 

111.3 

132.9 

122 . 9 

131.2 

71 . 6 

85.3 

109.7 

86.3 

134.0 

174.5 
181 . 8 

92 . 6 

151 . 9 

117 . 9 

72 . 2 

36.0 

52 . 8 

70.4 

62 . 6 

98.0 

112.4 
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TABLE 9. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO ICRF/DCRF 
MODELS (RMS units: nT) (con.) 





JULIAN 

RMS 

RMS 

RMS 

RMS 

PROJECT 

FLIGHT 

YEAR 

PAY 

X 

1 

z. 

£ 

C32-651 

5054 

1985 

321 

83.8 

129.6 

166.1 

145.2 

C32-651 

5055 

1985 

324 

50 . 6 

68 . 9 

74 . 4 

74 . 2 

C32-651 

4078 

1985 

327 

78 . 6 

133.7 

68 . 8 

68 . 6 

C32-651 

4079 

1985 

331 

78.6 

152.5 

82 . 6 

74.5 

C32-651 

4080 

1985 

334 

151.6 

247 . 3 

182 . 1 

16 6.7 

C32-651 

4081 

1985 

340 

89.5 

216.9 

119.7 

89 . 5 

C32-651 

4083 

1985 

343 

78.5 

96.0 

136.6 

143 . 6 

C32-651 

4084 

1985 

345 

64.8 

77.3 

75 . 3 

60.7 

C32-652 

1128 

1986 

27 

52.7 

53.2 

47 . 8 

53 . 1 

C32-652 

2019 

1986 

32 

62.2 

54.5 

95 . 4 

78 . 1 

C32-652 

5059A 

1986 

36 

49.6 

39.4 

88 . 6 

63 . 5 

C32-652 

5058 

1986 

40 

104.6 

42 . 9 

137 . 4 

67.5 

C32-652 

5057 

1986 

42 

84.3 

58.1 

154 . 9 

92 . 3 

C32-652 

5067 

1986 

47 

55.3 

61.7 

122 . 5 

102 . 3 

C32-652 

5068 

1986 

51 

65.3 

61.6 

127 . 9 

102.5 

C32-652 

5056 

1986 

53 

99.4 

52.0 

128 .4 

86.7 

C32-652 

5059B 

1986 

65 

63.7 

58.7 

116.9 

61 . 4 

C32-652 

5059C 

1986 

71 

54.5 

56.6 

70 . 1 

64 . 3 

C32-652 

5059D 

1986 

77 

44.2 

42.8 

58 . 7 

41 . 9 

C32-652 

5060 

1986 

79 

45.0 

71.0 

69 . 4 

49 . 8 

C32-652 

5061 

1986 

81 

37.0 

72.0 

94.0 

35 . 8 

C32-652 

5062 

1986 

85 

54.2 

70.0 

78 . 4 

54.2 

C32-652 

5063 

1986 

88 

100.8 

74.4 

150.5 

139.5 

C32-652 

5065 

1986 

95 

102.7 

75.0 

122 . 6 

122.2 

C32-652 

4085 

1986 

98 

57.0 

70.3 

103.1 

95 . 9 

C32-652 

3123 

1986 

101 

56.3 

54.8 

65 . 6 

65.2 

C32-652 

3124 

1986 

104 

51.9 

99.5 

82.4 

82 . 6 

C32-754 

3147 

1987 

261 

135.5 

116.5 

90 . 1 

98 . 6 

C32-751 

3142 

1986 

297 

75.0 

72 . 6 

71 . 5 

74.4 

C32-751 

3130 

1986 

327 

85.1 

52 . 9 

83.4 

77.5 

C32-751 

3131 

1986 

329 

107.6 

51.1 

81 . 1 

103.2 

C32-751 

3132 

1986 

334 

89.2 

50.9 

95.7 

85 . 8 

C32-751 

3133 

1986 

337 

77.2 

55.6 

71.5 

71 . 6 

C32-751 

3140 

1986 

343 

69.1 

72.8 

63.6 

60 . 8 

C32-751 

3134 

1986 

345 

55.3 

58.7 

66.1 

53 . 7 

C32-751 

3141 

1986 

347 

40.7 

48 . 9 

64 . 4 

51 . 4 

C32-751 

3135 

1986 

350 

66.6 

123.2 

143 . 7 

52.2 

C32-751 

3136 

1986 

353 

44.1 

70.3 

61 . 1 

39 . 8 

C32-751 

3138 

1986 

354 

53.4 

51.6 

55 . 9 

63.5 

C32-753 

1154 

1987 

108 

107.4 

101.8 

140.2 

138 . 1 

C32-753 

6019 

1987 

147 

100.2 

125.5 

130.4 

128.4 

C32-753 

1164 

1987 

181 

128.8 

120.9 

135.7 

136.7 

C32-753 

1160 

1987 

154 

90.1 

74 . 1 

117.5 

98.5 

C32-753 

1161 

1987 

156 

31.1 

76.9 

15.7 

21 . 4 

C32-754 

1159 

1987 

203 

69.9 

52.8 

71 . 6 

63 . 2 

C32-754 

3143 

1987 

206 

92.5 

59.7 

80.2 

94 . 0 

C32-754 

3144 

1987 

212 

74.6 

52.2 

78.7 

79.4 

C32-754 

4098 

1987 

216 

67.6 

135.7 

97.0 

65 . 1 

C32-754 

3145A 

1987 

231 

32.8 

200.2 

100.4 

36 . 8 

C32-552 

1133 

1985 

30 

88.6 

260.8 

138 . 9 

124.5 

C32-754 

5069 

1987 

233 

83.8 

181.5 

85.2 

88 . 6 
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TABLE 


PROJECT 

C32-754 

C32-754 

C32-754 

C32-754 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-851 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-852 

C32-952 

C32-952 

C32-952 

C32-952 

C32-553 

C32-651 

C32-651 

C32-352 

C32-853 

C32-453 

C32-453 

C32-951 

C32-951 

C32-951 

C32-951 

C32-951 

C32-951 

C32-951 

C32-951 

C32-951 


. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO IGRF/DGRF 
MODELS (RMS units: nT) (con.) 


FLIGHT 

5071 

3145 
4100 

3146 

3148 

3149 

3150 
3152 

3154 

3155 

3156 

3157 

3158 

3159 

3160 

3161 
1167 

4101 

4102 

4103 

4104 

4106 

4107 
4111 

4108 

3162 

3163 
3167 
4110 
4105 

3164 

1177 

1178 

1179 

1180 
3115 
3123A 
3124A 

5044 

4113 

4063 

4064 

4122 

4123 

4124 

4125 

4126 

4127 
4120 
4135 
5084 


YEAR 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1987 

1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 

1988 

1989 
1989 
1989 
1989 
1985 
1985 
1985 

1983 
1988 

1984 
1984 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 
1988 


JULIAN 

GM. 

251 

245 
251 

246 
291 
295 
299 
306 
312 
316 
320 
325 

29 
32 

335 

348 

20 

24 

27 

30 
34 
39 
41 
43 

49 

50 
56 
59 
63 
62 
77 
77 
79 
81 
83 

151 
292 

294 
69 

153 
150 

152 
283 
291 

295 
295 
303 
310 
297 
345 
324 


RMS 

X 

43.4 

63.3 

64 . 6 

116 . 9 

98 . 6 

61.9 

103.3 

130.9 

69.6 

76.5 

66.1 

81.5 

146.3 

93.8 

63.7 

78 . 8 

85.2 

42.1 

45.5 

51.7 

103.8 

78.9 

53.9 

84.9 

83.4 

71 . 9 

58.3 

78.4 

49.2 

71 . 8 

87.9 

118.1 

115.5 

111.8 

138.4 

122 . 6 

53.2 

56.3 

97 . 6 

80.3 

71.3 

55.7 

100.9 

96 . 6 

77 . 4 

130.2 

103 . 0 

57 . 1 

107.7 

71 . 1 

78.8 


RMS 

1 

54.7 
109.9 

70 . 3 

81 . 8 
10C .0 

104.6 

88.3 

97.7 
131.4 

60.7 
7 6.6 

159.6 

87 . 1 

78 . 8 
68 . 6 
45.3 

58.2 
32 . 6 

67 . 9 

59.1 

40.7 

70.3 

63.7 

66.2 

56.2 

54 . 9 

74 . 1 
131.6 

59.8 

72.3 
49.7 

98.2 
88 . 6 

79.2 
95.6 

142.2 

162.1 

188 . 8 

88.2 

89.1 

60.2 

53.8 
135.6 

67.8 
66.2 

82.5 
129.4 

77 . 4 
69.3 

39.9 

44 . 6 


RMS 


128.0 

78 . 5 
134 .2 

114.4 

103.0 

96.1 

104 . 1 

116 . 4 

77 .5 
68 . 1 

81 . 6 

112 . 9 

94 . 3 

110 . 4 

75 . 5 

78 . 0 

101 . 3 

84 . 9 

61.1 

68 . 3 
97 . 6 

80.6 

80.4 

65.7 

93.5 

82.7 

116 . 4 

105.8 

94.5 
67 . 4 

113.9 

139.4 

136.5 

119.8 

122 . 8 

113.2 

88.3 

86.5 

92.9 

114 . 8 

108.6 

87 . 6 

151.5 

106.5 

91.0 

75.5 

114 . 6 

82.7 

80.0 

79 . 4 

108.0 


RMS 

£ 

62.5 

59.7 

61.8 
90 . 3 

105 . 9 
80 . S 

111.4 

108.3 

69.0 
74 . 5 

66.7 
88 . 9 

100.8 

55.5 

73.9 

60.5 

71.7 

33.3 
47 . 8 
44 . 3 

57 . 6 

77.7 

40 . 2 

70 .2 
81 . 4 

74 .2 
67 . 6 

71.4 

56.6 

63.2 

59 . 9 

147.2 

147.5 

133.8 

149.6 

88.0 

63.0 

64 . 0 

67.7 
87 . 9 

78 . 8 

58.4 

160.6 

92 . 8 

80.2 

122 . 6 

92.5 

73.5 

91.0 

48.7 

81.1 
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TABLE 9. PROJECT MAGNET FLIGHT STATISTICS RELATIVE TO IGRF/DGRF 
MODELS .(RMS units: nT) (con.) 





JULIAN 

RMS 

RMS 

RMS 

RMS 

PROJECT 

FLIGHT 

YEAR 

DAY 

' X 

1 

z 

£ 

C32-253 

1079 

1982 

145 

59.4 

53.4 

80.5 

75 . 1 

C32-853 

5073 

1988 

117 

117.5 

93.1 

261 . 8 

263.4 

C32-853 

5081 

1988 

119 

73 . 9 

154 . 1 

190.2 

182.8 

C32-853 

5082 

1988 

122 

114.2 

110.7 

209 . 3 

182.8 

C32-853 

1168 

1988 

105 

45.9 

99.5 

66.4 

58.6 

C32-853 

1169 

1988 

174 

82.2 

57.7 

106.3 

58.0 

C32-853 

5076B 

1988 

131 

72.2 

98 . 1 

86 . 3 

84.0 

C32-853 

5083 

1988 

124 

133.5 

60.2 

257 . 0 

189.4 

C32-853 

4114 

1988 

164 

61.5 

76.7 

85 . 1 

7 4.5 

C32-853 

4115 

1988 

168 

108 . 9 

96.5 

99 . 6 

87 . 3 

C32-853 

4116 

1988 

171 

104 . 9 

66.0 

91 . 4 

88 . 8 

C32-951 

1171 

1988 

277 

74.8 

55.5 

79.4 

78 . 2 

C32-951 

3177 

1988 

307 

62.5 

79.5 

96 . 4 

87 . 1 

C32-951 

4119 

1988 

281 

70 . 3 

71.7 

75.5 

63.4 

C32-951 

4121 

1988 

299 

113.5 

117.9 

195.2 

215.4 

C32-951 

4128 

1988 

314 

73.2 

103.4 

85 . 3 

84.8 

C32-951 

4129 

1988 

318 

82 . 6 

81.7 

140 . 2 

135 . 7 

C32-951 

4130 

1988 

321 

121 . 6 

97.9 

107.2 

112 .2 

C32-951 

4131 

1988 

328 

73 . 1 

90.9 

141 . 5 

126.7 

C32-951 

4132 

1988 

331 

52.1 

89.6 

140 . 1 

109 . 6 

C32-951 

4133 

1988 

333 

70 . 6 

70.4 

66 . 9 

68 . 8 

C32-954 

1184 

1989 

207 

104.5 

194.1 

101.5 

109.3 

C32-954 

1186 

1989 

212 

106.9 

171 . 0 

99.2 

94 . 0 

C32-954 

1188 

1989 

218 

149.6 

111.6 

94.7 

83.0 

C32-954 

1189 

1989 

223 

99.8 

136.8 

114.4 

111 . 8 

C32-954 

1190 

1989 

228 

218.5 

160.4 

166.2 

118.5 

C32-954 

1196 

1989 

270 

99.7 

72 . 6 

79.6 

72 . 7 

C32-954 

1197 

1989 

229 

173.3 

112.3 

129 . 4 

8 6.6 

C32-954 

1198 

1989 

234 

266.8 

79.3 

102.3 

83.3 

C32-954 

1199 

1989 

244 

108.4 

72.5 

66.0 

68 . 9 

C32-954 

4139 

1989 

255 

102.5 

116.5 

93.2 

78 . 1 

C32-954 

1195 

1989 

267 

52.4 

119.1 

104 . 3 

89.4 

C32-954 

2057 

1989 

263 

81.3 

40.7 

107 . 8 

95.3 

C32-954 

2052 

1989 

247 

103.0 

83.9 

161.4 

173.7 

C32-954 

2055 

1989 

258 

77 . 9 

68.3 

88.2 

93 . 2 

C32-954 

2056 

1989 

261 

55.2 

83.2 

173.0 

131 . 6 

C32-954 

7012 

1989 

252 

99.8 

64.4 

93 . 9 

45 . 3 
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to the uniformity of satellite data, the orbit-by-orbit statist 
for MAGSAT and DE-2 data were taken to be the same as for each ent 
data set for each magnetic component. Note that data occupying ce 
for which there were fewer than 10 points were discarded due to 
presumed unreliability of their cell statistics. 

Now, the double summation expression over degree n and order m ^ 
B (r, 9, <t>) , B e (r, 0,4>) and B,(r,0,4>) can be converted to single summat 
expressions over the coefficient number / which ranges from 1 to - 
Then, we may write: 


1CS 

ire 

11s 

the 


ion 


B,(r,0,<|>) = 

1UO 

/ = l 

(34a) 

B e (r,6,<D) = 

168 

1 C,Q*{rM) 

l = I 

(34b) 

B,(r,0,<D) = 

168 

/*l 

(34c) 

B f (r , 0, <t>) = 

168 

I C,Q Fl (r,Q,ty) 
/ = 1 

( 34d) 


where the set of coefficients are the Gauss coefficients g Hm 

and h nn arbitrarily arranged so that: 

f l(n,m) = n(n + l)/2 + m 1 (35) 

C ‘ ~ \ Km* l(n,m) = n(n-\)/2 + m+90 J 


This ordering then requires: 




( n Y + 2 


(n + 1) 

— cosm<J>P7(0) ; l(n,m)- n(n + l)/2 + m 

Qrl( r ’ Q >® * ' 


\ r ) 

( o Y +2 ~~ 


(n + 1)1 

-5J sinm<}>P?(0) ; l(n,m) = n(n - l)/2 + m +90 




r 


Re 



Y+2 

/ 


rfPT(9) 

c0 sm *~dT '■ 

. rfPTO) 

mm *-dT '■ 


/(n,m) = n(« + 1)/2 + m 
l(n,m) = n(n - l)/2 + m +90 


(36a) 


(36b) 
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n + 2 


; l(n,m) = n(n + l)/2 + w 


(2^(r, 0, <|)) 



-m 



sin m4>P7(0)/ sin 0 
2 

cos /n 4>P7(0)/ sin 0 


V r ) 


l(n,m ) = n(n — 1 )/2 + m+ 90 f 


(36c) 


Qfi — 


1 

B f (r , 9, 0) 


{ B r (r,6, $)Q rl (r, 6, <p) + B e (r , 6, $)£?*(/•, 0, <J>) + B t (r, 0, 0)Q ti (r , 0, <D) } 


( 3 6 d ) 


These expressions are the most useful forms in which the spherical 
harmonic equations for the magnetic field components can be cast for 
a least-squares problem. 

Minimization of the chi-square function then requires that 


8X 2 


168 

I 


2 = 1 



(37) 


be a minimum, where the symbol 5 means variation. This in turn 
requires : 


dC, 


0 


j = 1 168 


(38) 


Therefore, we must have: 


dx? fto* &Xs 

dq ac, + ac, + ac, 


i. 


168 


(39) 


which is a nonlinear system of 168 equations for the 168 unknown 
coefficient set {C,}^. This system of equations is nonlinear since 

X 2 depends on B f which depends nonlinearly on the coefficients through 
the expression: 


= V B 2 r (r, 6, <j>) + B |(r , 9, + B J(r , 6, <J>) 

where, B,, B e , and B^ all depend linearly on the coefficients. 
Consequently, after noting that: 

ac, 


168 


1 c,l w n Q rl (r t A^ t )Qrj( r ,A^,) 


2 WrKQrMAA) 

i = 1 


(41a) 
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(41b) 


3 x 5 

dCj 

dC, 

9C, 


I c, I wMnMtevMM ' 

1=1 i =1 


168 

I 

(=1 1=1 


168 


X C, L "v. 
1=1 •=! 


i w Fl {Q rl {r„ QiA)QrM>’ % <*>.) + < 2«(0 0 i* < t>.)< 2 e/ r »> e » M + 


we have : 
168 


IQ2„ = % j = 1 ’-’ 168 


1 = 1 


where : 


and 


i M 1 

i^e*^.e^ i )e./ r '- 9 '>*' >+ l, WR{a,(ri,8 '^ ,)a ' (r, ’ e ''*' )+ 

i = 1 

8„ 0i)Qe>( r i> 0 *’ ♦() + ^M r i’ M h)Q*j( r i > e *’ W 

i >• 

*. * t W ri b ri Qrj<Ti >% 40 + .2, wAhQbM*’ e ‘> 40 + 

1 1=1 1-1 

I V 

I w ¥ b ¥ Q^(r it 0;, 40 + .i w F i b FiQFi( r i ’ 0 *’ 40 

be written in matrix form 


i * 1 


This system of 168 equations can 

cq =m " 

which has the inverse: 

C = O ' 1 * 


(410 

(41a) 

(42) 

(43) 

(44) 

as : 

(45) 

(46) 
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This is nsi the solution to the problem, however, since the right-hand 
side of this equation also depends on the unknown coefficients C t . 
That is, each element, %, of the vector ft depends on Q f (r,0,< », which 
depends on the unknown coefficients C, in a very nonlinear way. 

In order to solve this system of equations, we must iterate. If p 
is the iteration index, then we can let: 

c ( p ) = c, (p_,) + 5C < 

Then, in matrix form, we choose the following iteration scheme. 

C®> = C-'#’-" P = 1.2.. -P™, <43> 


The maximum numbe- -i iterations p„, is determined by requiring that: 

m I 491 

£5 C t < 1 nano tesla at p = 

i* i 


This condition must be tested after each iteration until it is 
satisfied • 

The rate of convergence of the algorithm depends strongly on the 
amount of noise (i.e., crustal influences, etc.) in the data. 
Aeromagnetic data have a great deal of crustal noise m them. Filtering 
the data to remove short wavelength 700 km) features from the data 
can improve the convergence rate by an order of magnitude. However, 
it has been shown that one-dimensional filters along the survey -rack 
leave short wavelength biases in the cross-track direction which 
adversely affect the final model. Consequently, no filtering (except 
for a very short wavelength 7 km) low-pass filter) was done on the 
aeromagnetic data. The number of iterations required for the model 

was p m „ = 14 . 

In order to implement the algorithm it is necessary to have an initial 

guess solution C (0) that is as close as possible to the actual solution. 
The a priori model coefficients used were the existing WC-85 model 
coefficients (Quinn, Kerridge, and Barraclough, 1986) . 

Notice that- in. this iteration scheme, the Q matrix, which has (168 
x 168) elements (as does its inverse Q~ l ) , needs to be computed only 
once since it does not depend on the coefficients C,. Note, too, that 
Q is a symmetric matrix so that only half of the elements in Q 
actually need to be computed. 
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SECTION 3. DISCUSSION 
3.0 Modeling Results 

An indication of the erratic nature of the geomagnetic field is 
provided by the wandering of the North and South Geomagnetic Poles. 
The pole movements since 1945 are illustrated in charts 5 and 6 which 
are based on the International Geomagnetic Reference Field ( IGRF ) 
models, WC-85 (revised) and WMM-90. The pole movements illustrate 
a poorly understood phenomena known as the geomagnetig. jerk whi :h 
occurred around 1970. The South Magnetic Pole movement in particular 
illustrates a sudden change in direction at about that time. T..ese 
jerks occur only a few times per century and are thought to be due 
to a sudden release of magnetic energy built up from the electromagnetic 
coupling between the top of the fluid core and the lower mantle, both 
of which have substantial electrical conductivities. The numerical 
pole positions at one _ year intervals for both poles are listed in 
table 10. 

It should be noted that roughly 90 percent of the Earth s magnet-*-*-' 
field is contained in the degree 1 spherical harmonic coefficients. 

£,°, g, 1 , and h\ . These three coefficients characterize the Earth's 
magnetic dipole field and form the basis of the geomagnetic coordinate 
system, which for the 1990 epoch is illustrated in chart 7. The axis 
of the geomagnetic coordinate system pierces the Earth' s surface at 
the Earth's magnetic dipoles, which are different from the dip poles. 
The location of the dipole poles is determined when the horizontal 
(H) component of the dipole field, computed from just the degree n = 1 
coefficients, is equal to zero. The dip poles on the other hand are 
determined when the horizontal ( H ) component of the field is computed 
using all 168 coefficients of the full-degree n = 12 model is equal to 
zero. For the WMM-90 model at 1990.0, the North magnetic dipole pole 
position is located at +79.35 degrees latitude and -71.10 degrees 
longitude, while the South magnetic dipole pole position is located 
at -79.35 degrees latitude and +108.86 degrees longitude. The 
displacement vector for the eccentric dipole for 1990.0 in the usual 
Earth-fixed spherical coordinate system (i.e., Z-axis is the rotation 
axis, X-axis points to the prime meridian and the X-axis is orthogonal 
to the other two, thereby creating a right-handed system) is 512 km 
radially outward from the Earth's center, with a colatitude 21.12 
degrees and a longitude of 145.70 degrees. 

A grid of main field and annual change values of the Earth's magnetic 
field derived from WMM-90 are tabulated in table 11 for seven basic 
magnetic field components ( X,Y,Z,H,F,D,l ) . Contours of five of these 
components (Z,H ,F ,D ,1 ) for the main field are illustrated in charts 
8 through 12. Contours of the annual change of these five components 
are illustrated in charts 13 through 17. These charts were plotted 
on a corrected Mercator projection. 
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CHART 5. NORTH MAGNETIC POLE MOVEMENT 
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CHART 6. SOUTH MAGNETIC POLE MOVEMENT 



TABLE 10. 


DIP POLE POSITION'S 


Year 

1 1945.000 

2 1946.000 

3 1947.000 

4 1948.000 

5 1949.000 

6 1950.000 

7 1951.000 

8 1952.000 

9 1953.000 

10 1954.000 

]_]_ 1955.000 

12 1956.000 

13 1957.000 

14 1958.000 

15 1959.000 

16 1960.000 

17 1961.000 

18 1962.000 

19 1963.000 

20 1964.000 

21 1965 . 000 

22 1966.000 

23 1967.000 

24 1968 . 000 

25 1969.000 

26 1970.000 

27 1971.000 

28 1972.000 

29 1973.000 

30 1974.000 

31 1975.000 

32 1976.000 

33 1977.000 

34 1978.000 

35 1979.000 

36 1980.000 

37 1981.000 

38 1982.000 

39 1983.000 

40 1984.000 

41 1985 .000 

42 1986.000 

43 1 9£T7 . 000 

44 1988.000 

45 1989.000 

46 1990.000 

47 1991.000 

48 1992.000 

49 1993.000 

50 1994.000 

51 1995.000 


North Pole 


Latitude 

(degrees) 

73 . 90 

74.05 
74 .20 

74 . 35 

74.50 

74 . 65 

74.75 

74.85 
74.95 

75.05 

75.20 

75.20 

75.25 

75.25 

75.30 

75.30 

75.35 
75.45 

75.50 
75.55 

75.60 

75.65 
75.70 

75.75 
75 .80 

75.90 
75 . 95 

76.00 

76 . 05 

76.10 
76 . 15 

76.30 
76.40 
76.55 

75.65 

76.90 

77.00 

77 . 10 

77.20 

77.30 

77 . 60 

77.75 

77.85 

78.00 

78 . 10 

78.25 

78.35 

78.50 

78 . 60 
78.70 

78 . 85 


Longitude 

(degrees) 

-100.20 
-100.35 
-100.45 
-100 . 60 
-100.75 
-100 . 85 
-100 . 95 
-101.10 
-101.20 
-101.25 
-101.45 
-101 . 35 
-101 . 30 
-101.20 
-101.15 
-101.05 
-101.10 
-101.15 
-101.25 
-101.25 
-101.35 
-101.25 
-101.20 
-101.10 
-101.05 
-101.00 
-100.90 
-100.80 
-100.70 
-100.60 
-100.65 
-100.85 
-101.05 
-101.25 
-101.45 
-101.70 
-101.90 
-102.10 
-102 . 30 
-102.55 
-102 . 60 
-102.80 
-103.00 
-103.25 
-103.40 
-103.70 
-103.85 
-104 . 15 
-104.35 
-104.55 
-104.80 


South Pole " 


Latitude 
(degrees ) 

-68 .15 
-68 .10 
-68 .05 
-68.00 
-67 . 94 
-67 .89 
- 67.75 
-67 . 62 
-67 .48 
-67 . 34 
-67.20 
-67 .10 
-67 .00 
-66 . 91 
-66 .81 
- 66.70 
-66 . 63 
-66.55 
-66.48 
-66.41 
-66.33 
-66.27 
- 66.21 
-66.15 
-66.09 
-66.03 
-65 . 96 
-65 . 90 
-65 . 84 
- 65.77 
-65.74 
-65 . 69 
-65.63 
-65.58 
-65.52 
-65.43 
-65 . 37 
-65 . 31 
-65.26 
-65.20 
-65.07 
-65.02 
-64 . 96 
-64 .91 
-64 .86 
-64 .80 
-64.75 
-64 . 69 
-64 . 64 
-64.58 
-64.53 


Longitude, 
(degrees ) 


144 . 

42 

144 . 

25 

144 . 

08 

143 . 

O 

143 . 

—1 1 
i 1 

143 . 

53 

143 . 

12 

142 

7 2 

142 . 

31 

141 . 

an 

> ''-J 

141 . 

,50 

141 . 

.24 

140 

.99 

140 

. 75 

140 

. 5 0 

14C 

O 1 

. C . J_ 

140 

.08 

139 

.93 

139 

.79 

139 

. 67 

139 

. 51 


139.43 

139.44 
139.41 
139 . 38 
139.40 
139.40 
139.43 
139.46 
139.46 
139.51 
139.51 
139.49 
139.49 
139.46 
139 . 35 
139.33 

139 . 32 

139.33 
139 . 31 
139.06 
139.00 
138.93 
138.87 
138.81 
138.74 
138 .69 
139.61 
138.56 
138.48 
138.43 
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Contours of these same five main field magnetic components plus grid 
variation in the north polar region are given in charts 18 through 
23, while contours of their secular variations are given in charts 
24 through 29. Similarly, for the south polar region, the main field 
contours are given in charts 30 through 35, while the corresponding 
secular variations are given in charts 36 through 41. These polar 
charts were plotted on a polar stereographic projection. Both the 
Mercator and polar stereographic charts were generated with respect 
to the 1984 World Geodetic System (WGS-84) ellipsoid. 

3 . 1 Final Comments 

The Polar Orbiting Geomagnetic Survey (POGS) satellite was launched 
in April of 1990, too late to be used in the 1990 epoch model. WMM-90, 
having been derived from data sets independent of POGS, will be a 
useful tool for evaluating POGS data and vice versa. Initial 
quantitative comparisons between WMM-90 and the POGS data indicate 
excellent agreement between the two. The POGS data will be used to 
fabricate the 1995 epoch model. Furthermore, if the satellite remains 
operational for its maximum expected lifetime of three years, it will 
for the first time be possible to generate a secular variation model 
to the same degree and order as the main field (i.e., N = A/ = 12) . 

Looking toward the end of this century and beyond, efforts have been 
made to secure data for modeling purposes via the Defense Meteorological 
Satellite Program (DMSP) platform. Efforts are underway to secure 
scalar data from a boom-mounted POGS-type magnetometer on the S-15 
DMSP satellite. This data will support the Epoch 2000 WMM . Further 
out, efforts are being made to secure full vector magnetic capability 
from DMSP Block 6 satellites that will operate during the first quarter 
of the next century. 
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u. S. NAVAL OCEANOGRAPHIC OFFICE 

CHART 7. GEOMAGNETIC COORDINATES 



TABLE 11. WMM-90 MAIN FIELD AND ANNUAL CHANGE GRID VALUES 
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CHART 9. VERTICAL COMPONENT (Z) 
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CHART 10. TOTAL INTENSITY (F) 
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CHART 15. TOTAL INTHNSI'IY (/ ) 
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CHART 29. GRID VARIATION (CV) 






(degrees) 
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uuoouuuuuouuuuuuuouuuuuuuuuuuuuuuuuuuouuuuouuuuuo 




SUBROUTINE GEOMAG (GEOMAGNETIC FIELD COMPUTATION) 


******************************************************************************** 


WMM-90 is a proposed Defense Mapping Agency (DMA) standard product. 
For information on the use and applicability of this product contact: 

DIRECTOR 

DEFENSE MAPPING AGENCY/HEADQUARTERS 
ATTN: CODE PR 
8613 LEE FUGHWAY 
FAIRFAX, VA 22031-2137 


******************************************************************************** 


GEOMAG PROGRAMMED BY: 

JOHN M. QUINN 7/19/90 

GEOPOTENTIAL DIVISION, CODE GGM 

U.S. NAVAL OCEANOGRAPHIC OFFICE (NAVOCEANO) 

STENNIS SPACE CENTER (SSC), MS 39522-5001 

PHONE: COM: (601) 688-4252 

AV: 485-4252 

FAX: (601) 688-5701 


PURPOSE: THIS ROUTINE COMPUTES THE DECLINATION (DEC), 
INCLINATION (DIP), TOTAL INTENSITY (TI), AND 
GRID VARIATION (GV - POLAR REGIONS USING A POLAR 
STEREOGRAPHIC PROJECTION ONLY) OF THE 
EARTH’S MAGNETIC FIELD IN GEODETIC COORDINATES 
FROM THE COEFFICIENTS OF THE CURRENT OFFICIAL 
DEPARTMENT OF DEFENSE (DOD) SPHERICAL HARMONIC WORLD 
MAGNETIC MODEL (WMM-90). THE WMM SERIES OF MODELS IS 
UPDATED EVERY 5 YEARS ON 1 JANUARY OF THOSE YEARS 
WHICH ARE DIVISIBLE BY 5 (I.E., 1980, 1985, 1990, ETC.) 

BY THE U.S. NAVAL OCEANOGRAPHIC OFFICE IN COOPERATION 
WITH THE BRITISH GEOLOGICAL SURVEY (BGS) AND 
IS BASED ON GEOMAGNETIC SURVEY MEASUREMENTS FROM 
AIRCRAFT, SATELLITE, AND GEOMAGNETIC OBSERVATORIES. 


Q ********************************************************************** ********** 
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MODEL- THE WMM SERES GEOMAGNETIC MODELS ARE COMPOSED 
OF TWO PARTS: THE MAIN FELD MODEL, WHICH IS 
VALID AT THE BASE EPOCH OF THE CURRENT MODEL, AND 
A SECULAR VARIATION MODEL, WHICH ACCOUNTS FOR SLOW 
TEMPORAL VARIATIONS IN THE MAIN GEOMAGNETIC FELD 
FROM THE BASE EPOCH TO A MAXIMUM OF 5 YEARS BEYOND 
THE BASE EPOCH. FOR EXAMPLE, THE BASE EPOCH OF 
THE WMM-90 MODEL IS 1990.0. THIS MODEL IS THEREFORE 
CONSIDERED VALID BETWEEN 1990.0 AND 1995.0. THE 
COMPUTED MAGNETIC PARAMETERS ARE REFERENCED TO THE 
WGS-84 ELLIPSOID. 


Q******************************************************************************** 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ACCURACY: IN OCEAN AREAS AT THE EARTH’S SURFACE OVER THE 
ENTIRE 5- YEAR LIFE OF A DEGREE AND ORDER 12 
SPHERICAL HARMONIC MODEL SUCH AS WMM-90, THE RMS 
DECLINATION ERROR IS ESTIMATED TO BE < .5 DEGREES, 

AND THE RMS INCLINATION ERROR IS ESTIMATED TO BE 
< .5 DEGREES. ALSO, THE RMS TOTAL INTENSITY ERROR 
IS ESTIMATED TO BE < 200 NANOTESLAS. 

THE ACCURACY AT ANY GIVEN TIME OF ALL FOUR 

GEOMAGNETIC PARAMETERS DEPENDS ON THE GEOMAGNETIC 
LATITUDE. THE ERRORS ARE LEAST AT THE EQUATOR AND 
GREATEST AT THE MAGNETIC POLES. 

IT IS VERY IMPORTANT TO NOTE THAT A DEGREE AND 

ORDER 12 MODEL, SUCH AS WMM-90, DESCRIBES ONLY 

THE LONG WAVELENGTH SPATIAL MAGNETIC FLUCTUATIONS 
DUE TO EARTH’S CORE. NOT INCLUDED IN THE WMM SERIES 
MODELS ARE INTERMEDIATE AND SHORT WAVELENGTH 
SPATIAL FLUCTUATIONS OF THE GEOMAGNETIC FIELD 
WHICH ORIGINATE IN THE EARTH’S MANTLE AND CRUST. 
CONSEQUENTLY, ISOLATED ANGULAR RESIDUALS AT VARIOUS 
POSITIONS ON THE SURFACE (PRIMARILY OVER LAND, IN 
CONTINENTAL MARGINS, AND OVER OCEANIC SEAMOUNTS, 
RIDGES, AND TRENCHES) OF SEVERAL DEGREES MAY BE 
EXPECTED. ALSO NOT INCLUDED IN THE MODEL ARE 
NONSECULAR TEMPORAL FLUCTUATIONS OF THE GEOMAGNETIC 
FIELD OF MAGNETOSPHERIC AND IONOSPHERIC ORIGIN. 

DURING MAGNETIC STORMS, TEMPORAL FLUCTUATIONS CAN 
CAUSE SUBSTANTIAL DEVIATIONS OF THE GEOMAGNETIC 
FIELD FROM MODEL VALUES. IN ARCTIC AND ANTARCTIC 
REGIONS, AS WELL AS IN EQUATORIAL REGIONS, DEVIATIONS 
FROM MODEL VALUES ARE BOTH FREQUENT AND PERSISTENT. 
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IF THE REQUIRED DECLINATION ACCURACY IS MORE 
STRINGENT THAN THE WMM SERIES OF MODELS PROVIDE, 

THE USER IS ADVISED TO REQUEST SPECIAL (REGIONAL OR 
LOCAL) SURVEYS BE PERFORMED AND MODELS PREPARED BY 
NAVOCEANO, WHICH OPERATES THE PROJECT MAGNET 
AIRCRAFT AND THE POLAR ORBITING GEOMAGNETIC SURVEY 
(POGS) SATELLITE. REQUESTS OF THIS NATURE SHOULD 
BE MADE THROUGH DMA AT THE ADDRESS ABOVE. 


******************************************************************************** 


USAGE: THIS ROUTINE IS BROKEN UP INTO TWO PARTS: 

A) AN INITIALIZATION MODULE, WHICH IS CALLED ONLY 
ONCE AT THE BEGINNING OF THE MAIN (CALLING) 

PROGRAM 

B) A PROCESSING MODULE, WHICH COMPUTES THE MAGNETIC 
FIELD PARAMETERS FOR EACH SPECIFIED GEODETIC 
POSITION (ALTITUDE, LATITUDE, LONGITUDE) AND TIME 

INITIALIZATION IS MADE VIA A SINGLE CALL TO THE MAIN 
ENTRY POINT (GEOMAG), WHILE SUBSEQUENT PROCESSING 
CALLS ARE MADE THROUGH THE SECOND ENTRY POINT (GEOMG1). 
ONE CALL TO THE PROCESSING MODULE IS REQUIRED FOR EACH 
POSITION AND TIME. 

THE VARIABLE MAXDEG IN THE INITIALIZATION CALL IS THE 
MAXIMUM DEGREE TO WHICH THE SPHERICAL HARMONIC MODEL 
IS TO BE COMPUTED. IT MUST BE SPECIFIED BY THE USER 
IN THE CALLING ROUTINE. NORMALLY IT IS 12 BUT IT MAY 
BE SET LESS THAN 12 TO INCREASE COMPUTATIONAL SPEED AT 
THE EXPENSE OF REDUCED ACCURACY. 

THE PC VERSION OF THIS SUBROUTINE MUST BE COMPILED 
WITH A FORTRAN 77 COMPATIBLE COMPILER SUCH AS THE 
MICROSOFT OPTIMIZING FORTRAN COMPILER VERSION 4.1 
OR LATER. 


****************************************************************** 


************** 


REFERENCES: 

JOHN M. QUINN, DAVID J. KERRIDGE, AND DAVID R. BARRACLOUGH, 
WORLD MAGNETIC CHARTS FOR 1985 - SPHERICAL HARMONIC 
MODELS OF THE GEOMAGNETIC FIELD AND ITS SECULAR 
VARIATION, GEOPHYS. J. R. ASTR. SOC. (1986) Vol. 87, 

PP. 1143-1157 
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DEFENSE MAPPING AGENCY TECHNICAL R EPORT TR 8350.2: 
DEPARTMENT OF DEFENSE WORLD GEODETIC SYSTEM 1984, 
SEPT. 30 (1987) 


JOSEPH C. CAIN, ET AL., A PROPOSED MODEL FOR THE 

INTERNATIONAL GEOMAGNETIC REFERENCE FIELD - 1965, 
I OF.OMAG. AND GEOELECL VOL. 19, NO. 4, PP. 335-355 
(1967) (SEE APPENDIX) 


ALFRED J. ZMUDA, WORLD MAGNETIC SURVEY 1957-1969, 
INTERNATIONAL ASSOCIATION OF GEOMAGNETISM AND 
AERONOMY (IAGA) BULLETIN #28, PP. 186-188 (1971) 


£******************* 


************************************************************* 


C 

C 

c PARAMETER DESCRIPTIONS: 


C 

C A 

C B 

C RE 

c SNORM 

c c 

C CD 

C TC 

c OTIME 

C OALT 

C OLAT 

C OLON 

C TIME 

C 

C ALT 

C GLAT 

C GLON 

C EPOCH 
C DTR 

C SP(M) 

C CP(M) 

C ST 

C CT 

C R 

C CA 

C SA 

C BR 

C BT 

C BP 

C P(N,M) 

C PP(N) 

C DP(N,M) 
C BX 

C BY 

C BZ 

C BH 

C DEC 

C 
C 


• SEMIMAJOR AXIS OF WGS-84 ELLIPSOID (KM) 

• SEMIMINOR AXIS OF WGS-84 ELLIPSOID (KM) 

■ MEAN RADIUS OF IAU-66 ELLIPSOID (KM) 

■ SCHMIDT NORMALIZATION FACTORS 

• GAUSS COEFFICIENTS OF MAIN GEOMAGNETIC MODEL (NT) 

- GAUSS COEFFICIENTS OF SECULAR GEOMAGNETIC MODEL (NT/YR) 

- TIME ADJUSTED GEOMAGNETIC GAUSS COEFFICIENTS (NT) 

- TIME ON PREVIOUS CALL TO GEOMAG (YRS) 

- GEODETIC ALTITUDE ON PREVIOUS CALL TO GEOMAG (YRS) 

- GEODETIC LATITUDE ON PREVIOUS CALL TO GEOMAG (DEG.) 

- GEODETIC LONGITUDE ON PREVIOUS CALL TO GEOMAG (DEG.) 

- COMPUTATION TIME (YRS) 

(E.G., 1 JULY 1985 = 1985.5G 

- GEODETIC ALTITUDE (KM 

- GEODETIC LATITUDE (DEG.) 

- GEODETIC LONGITUDE (DEG.) 

- BASE TIME OF GEOMAGNETIC MODEL (YRS) 

- DEGREE TO RADIAN CONVERSION 

- SINE OF (M*SPHERICAL COORD. LONGITUDE) 

- COSINE OF (M*SPHERICAL COORD. LONGITUDE) 

- SINE OF (SPHERICAL COORD. LATITUDE) 

- COSINE OF (SPHERICAL COORD. LATITUDE) 

- SPHERICAL COORDINATE RADIAL POSITION (KM) 

- COSINE OF SPHERICAL TO GEODETIC VECTOR ROTATION ANGLE 

- SINE OF SPHERICAL TO GEODETIC VECTOR ROTATION ANGLE 

- RADIAL COMPONENT OF GEOMAGNETIC FIELD (NT) 

- THETA COMPONENT OF GE . MAGNETIC FIELD (NT) 

- PHI COMPONENT OF GEOMAGNETIC FIELD (NT) 

- ASSOCIATED LEGENDRE POLYNOMIALS (UNNORMALIZED) 

- ASSOCIATED LEGENDRE POLYNOMIALS FOR M=1 (UNNORMALIZED) 

- THETA DERIVATIVE OF P(N,M) (UNNORMALIZED) 

- NORTH GEOMAGNETIC COMPONENT (NT) 

- EAST GEOMAGNETIC COMPONENT (NT) 

- VERTICALLY DOWN GEOMAGNETIC COMPONENT (NT) 

- HORIZONTAL GEOMAGNETIC COMPONENT (NT) 

- GEOMAGNETIC DECLINATION (DEG.) 

EAST=POSITIVE ANGLES 
WEST=NEGATIVE ANGLES 


(INPUT) 

(INPUT) 

(INPUT) 

(INPUT) 


(OUTPUT) 
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(OUTPUT) 



MAXDEG 

MOXORD 


- GEOMAGNETIC INCLINATION (DEG.) 

DOWN=POSITIVE ANGLES 
UP=NEGATIVE ANGLES 

- GEOMAGNETIC TOTAL INTENSITY (NT) 

- GEOMAGNETIC GRID VARIATION (DEG.) 

REFERENCED TO GRID NORTH 

GRID NORTH REFERENCED TO 0 MERIDIAN 
OF A POLAR STEREOGRAPHIC PROJECTION 
(ARCTIC/ ANT ARCTIC ONLY) 

- MAXIMUM DEGREE OF SPHERICAL HARMONIC MODEL 

- MAXIMUM ORDER OF SPHERICAL HARMONIC MODEL 


(OUTPUT) 

(OUTPUT) 


(INPUT) 


******************************************************************************* 


NOTE: THIS VERSION OF GEOMAG USES THE WMM-90 GEOMAGNETIC 
MODEL REFERENCED TO THE WGS-84 GRAVITY MODEL ELLIPSOID 


******************************************************************************* 


******************************************************************************** 


INITIALIZATION MODULE 


******************************************************************************** 


SUBROUTINE GEOMAG(MAXDEG) 


DIMENSION C(0: 1 2,0: 1 2) ,CD(0: 1 2,0: 1 2),TC(0: 1 2,0: 1 2) 
DIMENSION P(0:12,0:12),DP(0:12,0:12),SNORM(0:12,0:12) 
DIMENSION SP(0: 12),CP(0: 12),FN(0: 1 2),FM(0: 12),PP(0: 12) 
REAL K(0: 12,0: 12) 

EQUIVALENCE (SNORM.P) 


DATA EPOCH/1990.0/ 


203 



uuuuu 



DATA C/ 

0.0, 

-29780.5, 

-2134.3, 

1312.9, 

933.5, 


-208.3, 

* 


59.0, 

76.1, 

22.9, 

3.6, 

-3.3, 


1.3, 

* 


-1.3, 

5407.2, 

-1851.7, 

3062.2, 

-2244.7, 


784.9, 

* 


352.2, 

63.7, 

-62.1, 

2.3, 

9.5, 


-2.6, 

* 


-1.4, 

0.1, 

-2278.3, 

-384.3, 

1691.9, 


1246.8, 

* 


323.5, 

246.5, 

60.0, 

1.3, 

-1.2, 


-0.9, 

* 


4.5, 

-2.5, 

0.5, 

-284.9, 

291.7, 


-352.4, 

* 


808.6, 

-421.7, 

-110.8, 

-181.3, 

30.2, 


-11.7, 

* 


-10.7, 

-5.6, 

3.2, 

0.7, 

249.4, 


-232.7, 

* 


91.3, 

-296.5, 

139.2, 

-162.3, 

0.4, 


4.7, 

* 


-17.5, 

10.7, 

-3.6, 

0.2, 

0.4, 


40.8, 

* 


148.7, 

-154.6, 

-67.6, 

97.4, 

- 37.2, 


15.4, 

At 


7.9, 

2.2, 

-3.2, 

3.9, 

-1.1, 


-0.2, 

At 


-14.7, 

82.2, 

70.0, 

-56.2, 

-1.4, 


24.6, 

At 


-96.0, 

10.1, 

5.7, 

-1.4, 

3.2, 


0.3, 

At 


-1.1, 

-78.6, 

- 26.7, 

0.1, 

19.9, 


17.9, 

At 


-21.5, 

-6.8, 

1.9, 

3.0, 

6.3, 


1.7, 

* 


-0.3, 

0.9, 

9.7, 

-19.3, 

6.6, 


-20.1, 

* 


13.4, 

9.8, 

-19.0, 

-9.1, 

-7.0, 


0.8, 

At 


3.0, 

0.9, 

-0.6, 

-21.9, 

14.3, 


9.5, 

At 


-6.7, 

-6.4, 

9.1, 

8.9, 

-8.0, 


2.1, 

At 


-5.5, 

3.7, 

-1.1, 

0.8, 

2.6, 


1.2, 

At 


2.6, 

5.7, 

-4.0, 

-0.4, 

-1.7, 


3.8, 

At 


-0.8, 

-6.5, 

0.7, 

2.4, 

0.2, 


0.0, 

At 


1.0, 

-1.6, 

-2.2, 

1-1, 

- 0.7, 


-1.7, 

At 


-1.5, 

-1.3, 

-1.1, 

0.6, 

3.0, 


0.4, 

At 


0.7, 

0.7, 

1.3, 

-1.5, 

0.3, 


0.2, 

At 


-1.1, 

1.2, 

-0.2, 

-1.3, 

0.6, 


0.6, 

At 


0.2/ 








DATA CD/ 

0.0, 

16.0, 

-11.7, 2.1, 

-0.8, 

1.7, 

0.8, 

0.5, 

At 


0.0, 

0.0, 

0.0, 0.0, 

0.0, 

-13.8, 

9.3, 

3.7, 

At 


-7.6, 

1.0, 

0.0, 0.0, 

0.0, 

-1.1, 

0.0, 

0.0, 

At 


0.0, 

0.0, 

-12.8, -14.9, 

1.8, 

0.0, 

-7.4, 

0.0, 

At 


1.5, 

-0.9, 

0.0, 0.0, 

0.0, 

0.0, 

0.0, 

3.1, 

At 


0.8, 

-11.3, 

-5.8, 0.8, 

- 2.7, 

0.0, 

1.5, 

0.0, 

At 


0.0, 

0.0, 

0.0, 0.0, 

3.3, 

3.7, 

2.8, 

0.0, 

At 


-6.4, 

0.0, 

0.0, 2.7, 

-2.1, 

0.0, 

0.0, 

0.0, 

At 


0.0, 

0.0, 

-2.1, 1.2, 

1.2, 

0.6, 

3.0, 

0.0, 

At 


-1.0, 

0.0, 

0.0, 0.0, 

0.0, 

0.0, 

-0.6, 

-0.6, 

At 


0.0, 

-2.3, 

0.0, 0.0, 

0.0, 

0.0, 

1.0, 

0.0, 

At 


0.0, 

0.0, 

0.0, 0.6, 

0.8, 

0.0, 

0.0, 

0.0, 

At 


0.4, 

0.0, 

0.0, 0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

At 


0.4, 

-0.8, 

0.5, 0.3, 

0.5, 

0.0, 

-0.7, 

0.0, 

At 


57*0.0/ 
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INITIALIZE CONSTANTS 


IF (MAXDEG .GT. 12) MAXDEG=12 

MAXORD=MAXDEG 

PI=3. 14159265359 

DTR=PI/180.0 

SP(0)=0. 

CP(0)=1. 

P(0,0)=1. 

PP(0)=1. 

DP(0,0)=0. 

A=6378.137 

B=6356.7523142 

RE=6371.2 

A2=A**2 

B2=B**2 

C2=A2-B2 

A4=A2**2 

B4=B2**2 

C4=A4-B4 


CONVERT SCHMIDT NORMALIZED GAUSS COEFFICIENTS TO 
UNNORMALIZED 

SNORM(0,0)=1. 

DO 20 N=l,MAXORD 

SNORM(N,0)=SNORM(N- 1 ,0)*FLOAT(2*N- l)/FLOAT(N) 

J=2 

DO 10 M=0,N 

K(N,M)=FLOAT((N-l)**2-M**2)/FLOAT((2*N-l)*(2*N-3)) 

IF (M .GT. 0) THEN 

FLNMJ=FLOAT((N-M+ 1 )*J)/FLOAT(N+M) 
SNORM(N,M)=SNORM(N,M-l)*SQRT(FLNMJ) 

J=1 

C(M- 1 ,N)=SNORM(N,M)*C(M- 1 ,N) 

CD(M- 1 ,N)=SNORM(N,M)*CD(M- 1 ,N) 

END IF 

C(N,M)=SNORM(N.M)*C(N,M) 

CD(N,M)=SNORM(N,M)*CD(N,M) 

10 CONTINUE 

FN(N)=FLOAT(N+l) 

FM(N)=FLO AT (N) 

20 CONTINUE 
K(1,1)=0. 


OTEME=- 1000. 
OALT=-1000. 
OLAT=-1000. 
OLON=-1000. 
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return 


******************************************************************************** 


PROCESSING MODULE 


************** 


****************************************************************** 


ENTRY GEOMG 1 (ALT,GLAT,GLON,TIME,DEC,DIP,TI,GV) 


C 

C 


C 

C 


DT=TIME-EPOCH ___ 

IF (OTIME .LT. 0. .AND. (DT .LT. 0. .OR. DT .GT. 5.)) THEN 
PRTNT * * * 

’ WARNING - TIME EXTENDS BEYOND MODEL 5-YEAR LIFE SPAN 
CONTACT DMA FOR PRODUCT UPDATES:’ 


PRINT*, 
PRINT *, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
END IF 


DEFENSE MAPPING AGENCY’ 
ATTN: Code PRS’ 

8613 LEE HIGHWAY’ 

FAIRFAX, VA 22031-2137’ 
(703)285-9197, AUTOVON 356-9197’ 

EPOCH = ’.EPOCH 
TIME = ’.TIME 


RLON=GLON*DTR 

RLAT=GLAT*DTR 

SRLON=SIN(RLON) 

SRLAT=SIN(RLAT) 

CRLON=COS(RLON) 

CRLAT=COS(RLAT) 

SRLON2=SRLON**2 
S RL AT2=SRLAT* *2 
CRLON2=CRLON* * 2 
CRLAT2=CRLAT**2 
SP(l)=SRLON 
CP(l)=CRLON 

CONVERT FROM GEODETIC COORDS. TO SPHERICAL COORDS. 

IF (ALT .NE. OALT .OR. GLAT .NE. OLAT) THEN 
Q=SQRT(A2-C2*SRLAT2) 
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Q1=ALT*Q 

Q2=((Q1+A2)/(Q1+B2))**2 

CT=SRLAT/SQRT(Q2*CRLAT2+SRLAT2) 

ST=SQRT( 1 .0-CT**2) 

R2=ALT**2+2.0*Q1+(A4-C4*SRLAT2)/Q**2 

R=SQRT(R2) 

D=SQRT(A2*CRLAT2+B2*SRLAT2) 

CA=(ALT+D)/R 

SA=C2*CRLAT*SRLAT/(R*D) 

END IF 
C 
C 

IF (GLON .NE. OLON) THEN 

DO 40 M=2,MAXORD 

SP(M)=SP( l)*CP(M- 1 )+CP( 1 )*SP(M- 1 ) 

CP(M)=CP( 1)*CP(M- 1 )-SP( 1 )*SP(M- 1 ) 

40 CONTINUE 
END IF 
C 
C 

AOR=RE/R 

AR=AOR**2 

C 

C 

BR=0. 

BT=0. 

BP=0. 

BPP=0. 

C 

C 

DO 70 N=l,MAXORD 
AR=AR*AOR 
DO 60 M=0,N 
C 

COMPUTE UNNORMALIZED ASSOCIATED LEGENDRE POLYNOMIALS 
AND DERIVATIVES VIA RECURSION RELATIONS 
C 

IF (ALT .NE. OALT .OR. GLAT .NE. OLAT) THEN 
IF (N .EQ. M) THEN 
P(N,M)=ST*P(N- 1 ,M- 1) 

DP(N,M)=ST*DP(N- 1 ,M- 1 )+CT*P(N- 1 ,M- 1 ) 

GO TO 50 
END IF 

IF (N .EQ. 1 .AND. M .EQ. 0) THEN 
P(N,M)=CT*P(N- 1 ,M) 

DP(N,M)=CT*DP(N- 1 ,M)-ST*P(N- 1 ,M) 

GO TO 50 
END IF 

IF (N .GT. 1 .AND. N .NE. M) THEN 
IF (M .GT. N-2) P(N-2,M)=0.0 
IF (M .GT. N-2) DP(N-2,M)=0.0 


207 



on on 


c 


50 


P(N,M)=CT*P(N- 1 ,M)-K(N,M) *P(N-2,M) 
DP(N,M)=CT*DP(N-1,M)-ST*P(N-1,M)-K(N,M)*DP(N-2,M) 


END IF 
END IF 
CONTINUE 


TIME ADJUST THE GAUSS COEFFICIENTS 


IF (TIME .NE. OTIME) THEN 
TC(N,M)=C(N,M)+DT*CD(N,M) 

IF (M .NE. 0) THEN 

TC(M- 1 ,N)=C(M- 1 ,N)+DT*CD(M- 1 ,N) 

END IF 
END IF 

ACCUMULATE TERMS OF THE SPHERICAL HARMONIC EXPANSIONS 

PAR=AR*P(N,M) 

IF (M .EQ. 0) THEN 
TEMPI =TC(N,M)*CP(M) 

TEMP2=TC(N,M)*SP(M) 

ELSE 

TEMPI =TC(N,M)*CP(M)+TC(M-1,N)*SP(M) 
TEMP2=TC(N,M)*SP(M)-TC(M-1,N)*CP(M) 

END IF 

BT=BT- AR*TEMP 1 *DP(N,M) 

BP=BP+FM(M)*TEMP2*PAR 

BR=BR+FN(N)*TEMP1*PAR 

SPECIAL CASE: NORTH/SOUTH GEOGRAPHIC POLES 

IF (ST .EQ. 0.0 .AND. M .EQ. 1) THEN 
IF (N .EQ. 1) THEN 
PP(N)=PP(N-1) 

ELSE 

PP(N)=CT*PP(N- 1 )-K(N,M)*PP(N-2) 

END IF 

PARP=AR*PP(N) 

BPP=BPP+FM(M)*TEMP2*PARP 
END IF 


60 CONTINUE 
70 CONTINUE 


IF (ST .EQ. 0.0) THEN 

BP=BPP 

ELSE 

BP=BP/ST 

END IF 
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c 

ROTATE MAGNETIC VECTOR COMPONENTS FROM SPHERICAL TO 
GEODETIC COORDINATES 
C 

BX=-BT*CA-BR*SA 

BY=BP 

BZ=BT*SA-BR*CA 

C 

COMPUTE DECLINATION (DEC), INCLINATION (DIP), AND 
TOTAL INTENSITY (TI) 

C 

BH=SQRT(BX**2+BY**2) 

TI=SQRT(BH**2+BZ**2) 

DEC=ATAN2(BY,BX)/DTR 

DIP=ATAN2(BZ,BH)/DTR 

C 

COMPUTE MAGNETIC GRID VARIATION IF THE CURRENT 
GEODETIC POSITION IS IN THE ARCTIC OR ANTARCTIC 
(I.E. GLAT > +55 DEGREES OR GLAT < -55 DEGREES) 

C 

OTHERWISE, SET MAGNETIC GRID VARIATION TO -999.0 
C 

GV=-999.0 

IF (ABS(GLAT) .GE. 55.) THEN 

IF (GLAT .GT. 0. .AND. GLON .GE. 0.) GV=DEC-GLON 

IF (GLAT .GT. 0. .AND. GLON .LT. 0.) GV=DEC+ABS(GLON) 

IF (GLAT .LT. 0. .AND. GLON .GE. 0.) G V =DEC+GLON 
IF (GLAT .LT. 0. .AND. GLON .LT. 0.) GV=DEC-ABS(GLON) 

IF (GV .GT. +180.) GV=GV-360. 

IF (GV .LT. -180.) GV=GV+360. 

END IF 
C 
C 

OTIME=TIME 

OALT=ALT 

OLAT=GLAT 

OLON=GLON 

C 

C 

RETURN 

C 

C 

END 
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